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1 Introduction 



The method of QCD sum rules, developed more than twenty years ago by Shif- 
man, Vainshtein and Zakharov (SVZ)PJ has become a widely used working tool 
in hadron phenomenology.^ The advantages of this method are well known. In- 
stead of a model-dependent treatment in terms of constituent quarks, hadrons 
are represented by their interpolating quark currents taken at large virtuali- 
ties. The correlation function of these currents is introduced and treated in 
the framework of the operator product expansion (OPE), where the short- 
and long-distance quark-gluon interactions are separated. The former are cal- 
culated using QCD perturbation theory, whereas the latter are parametrized in 
terms of universal vacuum condensates or light-cone distribution amplitudes. 
The result of the QCD calculation is then matched, via dispersion relation, to 
a sum over hadronic states. The sum rule obtained in this way allows to cal- 
culate observable characteristics of the hadronic ground state. Inversely, the 
parameters of QCD such as quark masses and vacuum condensate densities 
can be extracted from sum rules which have experimentally known hadronic 
parts. What is also very important, the interactions of quark-gluon currents 
with QCD vacuum fields critically depend on the quantum numbers (spin- 
parity, flavor content) of these currents. Therefore, interpolating hadrons with 
currents, one is able ±o understand why the hadrons with different quantum 
numbers are not aliketl Numerous properties of hadrons with various flavor con- 
tent have been calculated by the sum rule method. The results are encouraging 
and in most cases reveal a remarkable agreement with the experimental data. 
Therefore, whenever one needs to determine an unknown hadronic parameter, 
the QCD sum rule prediction is among the reliable ones. 

However, the accuracy of this method is limited, on one hand, by the 
approximations in the OPE of the correlation functions and, on the other 
hand, by a very complicated and largely unknown structure of the hadronic 
dispersion integrals. The latter are usually approximated by employing the 
quark-hadron duality approximation. Consequently, the applicability of sum 
rules and the uncertainties of their predictions must be carefully assessed case 
by case. 

QCD sum rules are discussed in many reviews lTlIj (the list can be further 
continued) emphasizing various aspects of the method. In this review we try 
to present and update the subject in a more concise, i.e. close to encyclope- 
dic style, explaining the basics, presenting the most interesting applications, 
overviewing the recent developments and assessing the current and future po- 

b The mere fact that the original SVZ workQ has already got more than 2200 citations 
reflects the amount of papers employing QCD sum rules. 
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tential of sum rules. In several cases we compare the predictions of this method 
with the results of lattice QCD. Our purpose is to convince the reader that 
the analytical QCD sum rules are to a large extent complementary to the nu- 
merical lattice simulations of hadrons. It is a very special occasion to write 
a review on QCD sum rules in honor of B.L. Ioffe who contributed to this 
field with landmark results, considerably enlarging the spectrum of hadronic 
problems treated within this method. 

The content of the review is as follows. Section 2 is an elementary introduc- 
tion written for readers who have no experience in QCD sum rules and would 
like to grasp the basics of the method. In Section 3 we present an overview 
of the current status of many important applications and discuss the possible 
improvements and the perspectives for new investigations. In Section 4 we 
consider the light-cone version of QCD sum rules, explaining the idea, outlin- 
ing the derivation and presenting the main applications. Section 5 contains a 
summary. 

2 Understanding SVZ sum rules 

2. 1 Correlation function of quark currents 
The QCD Lagrangian 

£qcd = -\Gl v G^ v + Mi P ~ rn q )% , (1) 

where G^ v is the gluon field-strength tensor and ip q are the quark fields with 
different flavors: q = u,d, s, c, .., is discussed in detail in many chapters of this 
book. It is our common belief that this Lagrangian governs all properties of 
hadrons and hadronic processes. However, a direct use of Eq. (Q) and of the 
corresponding Feynman rules is possible only within the limited framework 
of perturbation theory. At least some of the quarks or gluons in a hadronic 
process have to be highly virtual. This condition guarantees the smallness of 
the corresponding effective quark-gluon coupling a s — 3s/47r and, thereby, a 
legitimacy of the perturbative expansion. Usually, high virtuality is achieved 
in a scattering of hadrons at large momentum transfer. However, even for 
these specially configured hard scattering processes, a perturbative calculation 
of quark-gluon Feynman diagrams is not sufficient, because the quarks par- 
ticipating in the hard scattering are confined inside hadrons. Hence, one has 
to combine the perturbative QCD result with certain wave functions or mo- 
mentum distributions of quarks in hadrons. To calculate these characteristics, 
one needs to know the QCD dynamics at distances of order of the hadron 
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Figure 1: Quark-antiquark creation and annihilation by the virtual photon in the electron- 
electron scattering. 

size: Rhadr ~ 1/Aqcdj the scale at which perturbation theory in a s is not 
applicable. 

In order to avoid long-distance problems, one could consider processes 
with no initial and final hadrons, and with all quarks propagating at short 
distances (during short times). Such configurations are not that hypothetical 
as it may seem. They are realized in nature, when the quark-antiquark pair 
is produced and absorbed by an external source, e.g. by a virtual photon in 
the electron-electron elastic scattering, as shown in Fig. [I]. The propagation of 
the intermediate quark-antiquark pair in this process adds a very small (order 
of a em — e 2 /4ir ) quantum correction to the cross section of e~e~ — * e~e~. 
Nevertheless, taken separately, the amplitude of the quark-pair creation and 
annihilation is an extremely useful object from the QCD point of view. The 
formal expression for this amplitude can be written as 

n^(q) = if d 4 x e iq - x (0\ T{j^x)j v (0)} |0) = (q^ - q 2 g^)Tl(q 2 ) , (2) 

where q is the four-momentum of the virtual photon with q 2 <0, = ip'y^ip is 
the colorless quark current with a given flavor ip = u,d,s,c, .. (for simplicity, 
we have omitted the electromagnetic coupling from this definition). In the 
amplitude (0) , the initial and final states contain no hadrons and are therefore 
identified with the vacuum state of QCD. The Lorentz-structure of the r.h.s.of 
Eq. (||) is dictated by the conservation of the electromagnetic current: d^j^ = 
0, so that the single invariant amplitude II(g 2 ) encodes all dynamical effects. 
The amplitude II^,, represents an important example of a two-point corre- 
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lotion function (correlator) of quark currents. If the four-momentum squared 
transfered to the quarks is large, Q 2 = —q 2 ^S> Aq CD , this correlation function 
turns into the genuine short-distance object we are looking for. This is simply 
because the integral in Eq. (||) is dominated by small spatial distances and 
time intervals: 

\x\ ~ X ~ l/VO* < Rhadr ■ (3) 

This condition can be directly inferred from a general analysis of Eq. (|J) in the 
case of massless quarks, an approximation justified for the light u, d, s quarks. 
After contraction of the Lorentz indices, the vacuum average in Eq. (|^) can 
only depend on the space-time interval x 2 ~ Xq — x 2 : 

(0 | T{Ux)nO)} I 0) = J dr e^ 2 f(r) , (4) 

where the Fourier transform of this functional dependence is introduced. Using 
the representation (^) and shifting the variable x, one obtains from Eq. 

3q 2 Il(q 2 ) =-i [dr [ d 4 x e lTX \ lQ2/iT f{r) . (5) 



The integrand on the r.h.s. of this equation is suppressed if at least one of the 
exponential functions rapidly oscillates. Therefore, dominant contributions to 
II(<7 2 ) stem from the regions where both r ~ 1/x 2 and r ~ Q 2 . To fulfill these 
two conditions simultaneously, one has to demand 

x 2 ~ l/Q 2 , (6) 

so that, at Q 2 — > oo the quarks propagate near the light-cone, x 2 ~ . This 
condition is necessary but not yet sufficient for the short-distance dominance. 
To demonstrate the latter it is convenient to choose the reference frame go = 0, 
so that q 2 = Q 2 , and the exponent in Eq. (^|) is simply equal to exp(— iq ■ x). 
Again, to avoid a fast oscillating integrand one needs 



\x\~l/VQ 2 , (7) 

which, together with Eq. (||), yields Eq. (^). Hence, at large Q 2 the quarks in 
the diagram in Fig. ^ propagate predominantly at short distances and during 
short time intervals. Due to the asymptotic freedom of QCD, the quark-gluon 
interactions are then suppressed. Therefore, as a first approximation, one may 
calculate the correlation function (^) representing virtual quarks by the free- 
quark propagators inferred directly from the Lagrangian (0). In the case of 
heavy quark currents (ip = c,b), the situation is even simpler because the 
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quark mass m c b 3> Aqcd introduces an intrinsic large energy scale. One 
has an asymptotically free quark-antiquark fluctuation already at small q 2 -C 
4 m 2 b . In this case, the characteristic distances in the correlation function are 
determined by the inverse heavy quark mass |af| ~ x ~ 1/(2 m C) &). 

2.2 Summing up hadrons: the unitarity relation 

Before turning to the actual calculation of II M „, let us discuss how this object 
is related to physically observed hadrons. Note that the invariant amplitude 
II((7 2 ) is an analytic function of q 2 defined at both negative (spacelike) and 
positive (timelike) values of q 2 and, formally, even at complex values of this 
variable. At positive q 2 the underlying electromagnetic process is the cross- 
channel of the electron-electron scattering, i.e. the annihilation e + e~ — > e + e~ 
with the total cm. energy E e + + E e - — \fq 2 ~. If q 2 is shifted from large 
negative to positive values, the average distance between the points and x 
in the quark amplitude (||) grows. The long-distance quark-gluon interactions 
become important and, eventually, the quarks form hadrons. In particular, a 
quark-antiquark pair created by the current j M with the spin-parity J p = 1~ 
materializes as a neutral vector meson. For the currents v/y^u and dj^d, one 
has to respect the isospin symmetry, because the mass difference md — m u ~ 
O(MeV) is much smaller than the QCD scale Aqcd- The ground-state vector 
mesons with the isospin 1 = 1 and 1 = are p and w with the quark content 
(uu — dd)/V2 and (uu + dd)/V2, respectively. The lightest vector meson 
created by the current sj^s is 4>. For the heavy quark currents C7 M c and 
&7 M 6, the ground states are J/tjj and T, respectively. Physically, vector mesons 
are observed in a form of resonances in e + e~ annihilation at energies \fq 2 = 




vciy , (V = p,oj,(f),J/ip,T,...). Not only the ground-state, but also excited 
vector mesons and a continuum of two- and many-body hadron states with 
the quantum numbers of V contribute to iI M1/ . 

A rigorous way to quantify a very complicated hadronic content of at 
q 2 > is provided by the unitarity relation obtained by inserting a complete 
set of intermediate hadronic states in Eq. (0) : 



where the summation goes over all possible hadronic states | n) created by 
the quark current including sums over polarizations, and dr n denotes the 
integration over the phase space volume of these states. 




2ImIV(< ? ) = ^(0 1 j M \n){n\j v |0) dr n (27r) 4 ^% - p n ) , 



(8) 



n 
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The one-particle, vector meson contribution to the hadronic sum (|8J) is 

-Im UY v (q 2 ) = (q^ - m v g^)fy8{q 2 - my) , (9) 

where the total decay width of V is neglected for simplicity, and the decay 
constant fy is defined by the matrix element of the current between the 
vacuum and the vector-meson states: 

(V(q)\j^\0) = f v m v e<y>, (10) 

e£i being the polarization vector of V (e v ■ q = 0). Note that fy is a typical 
hadronic parameter determined by the long-distance dynamics. 

The contributions of continuum hadronic states to the unitarity relation 
(^) are more involved. Each individual state \n) yields a continuous imaginary 
part at q 2 > m„, m n being the sum of hadron masses in this state. Moreover, 
the corresponding hadronic matrix elements for multiparticle states (n \ | 0) 
are not just constants but depend on q 2 . 

For convenience, we single out the ground-state vector-meson contribution 
on the r.h.s. of the unitarity relation (^) and introduce a compact notation 
for the rest of contributions including excited vector mesons and continuum 
states: 

-Im U(q 2 ) = f 2 S(q 2 - m 2 v ) + p h (q 2 )9(q 2 - S h ) , (11) 

where Sq is the threshold of the lowest continuum state. Notice that in the 
light quark channels this threshold, set by two- and three-pion states, is lower 
than my. In the heavy quark channels the pattern is different. There are 
several heavy quarkonium resonances below the threshold of the heavy fla- 
vored meson pair production and, therefore, Eq. ([ll]) has to be slightly modi- 
fied: Im II V (q 2 ) — > ^2 V Im H v (q 2 ) including the sum over all below-threshold 
quarkonium states. 

2.3 Deriving the dispersion relation 

From the above discussion we learned that the correlation function (||) is an 
object of dual nature. At large negative q 2 it represents a short-distance quark- 
antiquark fluctuation and can be treated in perturbative QCD, whereas at 
positive q 2 it has a decomposition in terms of hadronic observables. The next 
step is to derive a dispersion relation linking Tl(q 2 ) at an arbitrary point q 2 < 
to the hadronic sum (ph. For that, one employs the Cauchy formula for the 
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Figure 2: The contour in the plane of the complex variable q 2 = z. The open point indicates 
the q 2 < reference point of the QCD calculation. Positions of hadronic thresholds at q 2 > 

are indicated by crosses. 



analytic function Tl(q 2 ) 1 choosing the contour shown in Fig. |[ 

,9, i / , n(z) i / , n(z) 



2ni J z — q 2 2ni J z — q 2 

C ~ \z\=R 



27ri J z — q l 



The radius i? of the circle in this contour can be put to infinity. This simplifies 
the r.h.s. of Eq. ( p"2| ) considerably, because, if the correlation function vanishes 
sufficiently fast at |g 2 | ~ R — > oo, (if lim| g 2|^ oo n(g 2 ) ~ l/|g 2 | e , with any 
e > 0) the integral over the circle tends to zero. Below we shall discuss a 
necessary modification of Eq. ( p"2| ) if II(q 2 ) does not vanish. The second integral 
on the r.h.s. of Eq. ( p^ ) can be replaced by an integral over the imaginary 
part of U(q 2 ). One makes use of the fact that n(<? 2 ) is real at q 2 < t m i n = 
minim 2 /, Sp}. Hence, according to the Schwartz reflection principle: H(q 2 + 
ie) — n(q 2 — ie) = 2i ImIT(q 2 ) at q 2 > t m i n . After this replacement, we obtain 
the dispersion relation: 

oo 

ntf)= 1 - [ds^L. (is) 

7r J s — q A — ie 
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The infinitesimal — ie will not be shown explicitly hereafter. As we shall see 
in the next subsection, the correlation function (j^) is ultraviolet divergent. 
Consequently, the imaginary part Imil(s) does not vanish at s — > oo and the 
dispersion integral diverges. A standard way to cure this problem is to 
subtract from H(q 2 ) first few terms of its Taylor expansion at q 2 = 0. For the 
correlation function (||) one subtraction is sufficient: 

n( g 2 ) = n(< 7 2 )-n(o). (14) 

The dispersion relation (O) is modified in the following way: 



oo 

tt; 2n Q 2 f , ImII(s) , ir . 

n <r = — / ds - ( — -k ■ 15 

7T J s(s - q 1 ) 

tmin 

Using the hadronic representation (|Tl|), one finally obtains 

oo 
s h 

Notice that in our case 11(0) = 0, due to the gauge invariance of the electro- 
magnetic interaction. Nevertheless, we retain 11(0) in Eq. @, having in mind 
a generic case, where a subtraction constant or a finite polynomial in q 2 appear 
in the resulting dispersion relation. 

The dispersion relations similar to Eq. ( [l6| ) are central objects of our re- 
view. With the correlation functions calculated in QCD in a certain approx- 
imation, these relations establish sum rules, i.e. nontrivial constraints on the 
sums over hadronic parameters. 



2.4 Applying the Borel transformation 

The sum rules in the form (|l^) are not yet very useful, e.g. for estimating 
the parameters of the lowest-lying hadronic state. They are in general plagued 
by the presence of unknown subtraction terms. More importantly, little is 
known about the spectral function p h (sl of excited and continuum states. The 
situation can be substantially improved!^ if one applies to both sides of Eq. (16) 
the Borel transformation 

/ 2\(n+l) / J \" 

H(M 2 ) = B M ,U( q 2 ) = lim n(q 2 ). (17) 



— q 



-q 2 /n=M 2 



dq 
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Two important examples are 



n ( n 2 \* n n ( 1 ^ 1 exp(-m 2 /Af 2 ) 

at A: > 0. Jjpinsformations of more complicated functions can be found in the 
literature ffl Applying Eqs. (17) and (18) to Eq. (^6|), a more convenient form 
of the sum rule is obtained: 

/>oo 

n(M 2 ) = /2 e -</^ 2 + / p h («)e-/ Ma . (19) 



The Borel transformation removes subtraction terms in the dispersion relation 
and exponentially suppresses the contributions from excited resonances and 
continuum states heavier than V. In the case of the heavy quark-antiquark 
currents, instead of the Borel transformation, it is more useful to apply a 
simpler procedure of n-times differentiation of Eq. (|l^) at q 2 = q 2 < : 

i <** m os, _ p v , r,„ p h (s) 



One gains a power suppression of heavier states and again removes the sub- 
traction terms. 



2.5 Calculating the correlation function in QCD: the perturbative part 

We now turn to the next important step in the sum rule derivation and describe 
how the QCD calculation of the correlation function Q) is done. As explained 
in subsection 2.1, at very large Q 2 = —q 2 the function ILj„ can be approxi- 
mated by the free-quark loop diagram shown in Fig. |^a. In a generic case of a 
single quark flavor with the mass m, one has to contract all quark fields in (^|) 
considering the quark propagators in the free-quark approximation: 

S»(x,y) = -z(0|T{^(x)?(y)}|0) = & [ e -*-(x-v)^±IL , (21) 

J (27r) 4 — m z 

where the quark color indices are explicitly shown. After integrating over x, 
shifting to D ^= 4 dimensions and taking traces, we obtain 



q 



n(°V) = - L[J?P ( 2 - D ^ - W - + ° m2 . (22) 



(D - 1) J J (2tt) d (p2 + q 2 v{ i _ v - } _ m 2f 
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(a) 




(b) (c) (d) 

Figure 3: Diagrams determining the correlation function (Q): the free-quark loop (a), the 
perturbative QCD corrections (b,c,d). Solid lines denote quarks, dashed lines gluons, wavy 

line external currents. 



One may now proceed, performing the momentum integration in D dimensions. 
However, at this point we prefer to modify the standard procedure, and apply 
the Borel transformation before the integration in p is done. The divergence 
disappears and the limit D — > 4 can safely be restored. After the Wick rotation 
to the Euclidean space p 2 —> —p 2 = z and the angular integration in the four- 
dimensional integral, one obtains, in the massless quark case, 

oo 

B M2 (q 2 rt°Hq 2 )) = J dzf^(z), (23) 


where 

f{0} ^ = ^ Z j^){- 1 + ^^)^^ (24) 
o 

The numerical result for f^\z) is shown in Fig. [|. We see that the average (z) 
characterizing the quark virtuality in the loop diagram, is of the order of M 2 , 
and that the region of small z, e.g., z < 0.1 M 2 is strongly suppressed. Thus, 
at sufficiently large M 2 , the quarks are predominantly far off-shell. In the case 
of heavy c, b quarks it is more useful to differentiate Eq. ( B2n n times over q 2 
at q 2 = 0, revealing that the dominant contribution stems from \p 2 \ ~ m 2 b /n. 
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0.01 



Figure 4: The distribution of the quark virtuality in the loop diagram after the Borel trans- 
formation at M 2 = 1 GcV 2 (solid line) and M 2 = 10 GeV 2 (dashed line). 



Thus, an average high virtuality of heavy quarks is guaranteed, if n is not very 
large. 

The most convenient final expression of U(°\q 2 ) is in the form of the 
dispersion integral: 

tf°Hf) = £j*> ^$ (25) 



s(s — q 2 ) 



with the imaginary part 



Imn (0) (s) = — u(3 - v 2 )9(s - Am 2 ) . 
8n 



(26) 



where v = y/l — Am 2 / s. One should be careful in interpreting the imaginary 
parts of quark-loop diagrams. In QCD quarks are confined (the full quark 
propagators have no poles). Hence, the imaginary part (|2q ) is a purely math- 
ematical object. The free-quark approximation H^^q 2 ) is especially simple in 
the light-quark case m 2 Q 2 , yielding 



n(°V) 



ds _ 1 ln Q 2 | fm 2 
At: 2 J s(s-q 2 ) An 2 Am 2 \Q 2 

4m 2 



(27) 



where the 0(m /Q ) correction is numerically important only in the s-quark 
case. 
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To improve the free-quark approximation, one has to calculate the 0(a s ) 
perturbative correction corresponding to the diagrams in Figs. [|b,c,d. Since we 
already convinced ourselves that the average quark virtualities in the loop are 
of 0(M 2 ), it is conceivable to use QCD perturbation theory for these diagrams 
taking the quark-gluon coupling a s at the scale M. The calculation of two- 
loop diagrams is technically quite involved but, fortunately, the result can he 
directly taken from QED, employing the Schwinger interpolation formula Ej 
for the 0(a em ) radiative correction to the electron polarization operator. One 
has simply to replace a em — ► a s Ci? {Cp — 4/3). After adding the 0(a s ) 
correction, the perturbative part of the correlation function, ]j(P ert ') (q 2 ), is 
given by the dispersion relation fl2q) with the imaginary part 



Imfl 



(pert) 



Imn (0) (s)<! l + a s C F 



7T 

2^ 



3 / 7T 



3 

4tt 



(28) 



In the case of u,d,s quark currents, the 0(a s ) correction in ( B8[ ) reduces to 
a s /n and is numerically small. Virtual gluon exchanges are potentially impor- 
tant for the correlation functions of heavy quark-antiquark currents since the 
a s /v term in ( |28| ) becomes anomalously large at v <C 1, i.e., at s close to the 
threshold 4m 2 . The a s /v terms can be traced back to Coulomb- type interac- 
tions between quarks and antiquarks. In the nonrelativistic approximation, it 
is possibles to sum up all (a s /v) n terms in the correlation function, taking 
into account not only the one-gluon exchange but the whole ladder of such 
exchanges. This summation is usually applied to the correlation functions of 
&7 M 6 currents related to T resonances | 

A careful reader may have noticed that in all perturbative diagrams con- 
sidered here the regions of small quark and gluon virtualities are automatically 
included, e.g., the integration in ( |2~i| ) spreads over small quark virtualities. The 
QCD perturbation theory is invalid in this region and the free propagators can- 
not be used. Nevertheless, one may still argue that no large numerical error 
is being introduced, as long as M 2 is kept large. For instance, in the integral 
(E3T) the region with z < fi 2 contributes with a suppression of 0(/j> 4 /M 4 ). 



2. 6 Vacuum condensates and operator product expansion 

The fact that the perturbative part of II^,, has been reliably calculated does 
not yet imply that all important contributions to the correlation function have 

c At 0(a s ) accuracy one should take care of a proper definition of the heavy quark mass 
which is a scale-dependent parameter in perturbative QCD. In Eq. (GH) the so called "pole" 
mass is used. We do not discuss this important issue here, and refer to the chapters by 
Uraltsev and Chetyrkin in this book. 



14 



been taken into account. The complete calculation El has to include the effects 
due to the fields of soft gluons and quarks populating the QCD vacuum. The 
problem of vacuum fields in QCD with its many interesting aspects is discussed 
elsewhere in this book. We only mention that vacuum fluctuations in QCD are 
due to the complicated nonlinear nature of the Lagrangian ([l]) . The ultimate 
solution of QCD equations of motion and the resulting complete picture of the 
vacuum fields are unknown. Various nonperturbative approaches (instanton 
models, lattice simulation of QCD, etc.) indicate that these fields fluctuate 
with typical long-distance scales A vac ~ Aqcd- It is clear that the quark- 
antiquark pair created by the external current at one point and absorbed at 
another point interacts with the vacuum fields. This interaction is beyond 
QCD perturbation theory and has to be taken into account separately. 

A practical way to calculate the vacuum-field contributions to the cor- 
relation functions relies on the following qualitative arguments. At large 
Q 2 ^> Aq C£) , the average distance between the points of the quark-antiquark 
emission and absorption is essentially smaller than the characteristic scale of 
the vacuum fluctuations. Therefore, propagating in the QCD vacuum, the 
quark-antiquark pair acts as a short-distance probe of long-distance fields and 
perceives static, averaged characteristics of these fields. At the same time, 
the emission of quarks and antiquarks does not significantly disturb the vac- 
uum state. Hence, in the first approximation, quarks with large momenta 
~ y/Q 2 scatter over external static fields composed of soft vacuum gluons and 
quarks. The corresponding diagrams are shown in Fig. |5| In the case of light 
quarks, there are several important effects: the vacuum gluons are emitted 
and absorbed by virtual quarks (Figs. ||a,b,c), the quarks and antiquarks are 
interchanged with their vacuum counterparts (Fig. ||d) and, finally, a combined 
quark-gluon interaction takes place (Figs. |^e,f). For the heavy quarks only the 
interactions with the vacuum gluons are important. 

A quantitative framework which follows this picture and incorporates both 
short- and long-distance contributions was developed™ in a form of a general- 
ized Wilson OPE. To apply this method to the correlation function (|[), one 
has to expand the product of two currents in a series of local operators: 

i Jd 4 x e*- ^(1)7^(30,^(0)7,^(0)} 

= (q^-q 2 g^)^2C d (q 2 )O d , (29) 

d 

so that 

n(g 2 ) = £c- d (<z 2 )<0|O d |0) . (30) 

d 
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(a) 



(b) 



(c) 




(d) 



(e) 



(f) 



Figure 5: Diagrams corresponding to the gluon (a,b,c), quark (d), quark-gluon (e) and 
four-quark (f) condensate contributions to the correlation function (0). 



In this expansion, the operators are ordered according to their dimension d. 
The lowest-dimension operator with d = is the unit operator associated 
with the perturbative contribution: Co(s!) = W ert (q 2 ), (0 | O | 0} = 1. The 
QCD vacuum fields are represented in ( |30| ) in a form of the so called vacuum 
condensates, the vacuum expectation values of the d ^ operators, composed 
of quark and gluon fields, ip, tj) and G". The contributions of high-dimensional 
condensates corresponding to the diagrams with multiple insertions of vacuum 
gluons and quarks, are suppressed by large powers of A 2 ac /Q 2 . Therefore, even 
at intermediate Q 2 ~ 1 GeV 2 , the expansion (^0|) can safely be truncated after 
a few terms. 

One may interpret the OPE (|3(]) in the following way. The product of 
quark currents acts as a quasi- local "antenna" having a small size 0(\j \J Q 2 ) 
and probing the static vacuum fields. This interaction depends on the prop- 
erties of the "antenna" , i.e. on the quantum numbers and flavor content of 
the quark currents. In Eq. (|30|), this dependence is accumulated in the Wil- 
son coefficients Cd(q 2 ) receiving dominant contributions from the regions of 
short distances (large momenta). In addition, the current- vacuum interaction 
is determined by the long-distance dynamics represented in Eq. ( |30| ) by the 
universal condensates (0 | Od | 0) with d ^ 0, which are independent of the 
properties of the quark currents. The separation of distances is a key point in 
the OPE j3C|). Representing II (q 2 ) in this form, one introduces a certain scale 
ji which separates the regions of short and long distances. The interactions at 
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momenta p 2 > jj 2 are included in the coefficients Cd(q 2 ), while the effects at 
p 2 < fj, 2 are absorbed into the vacuum condensates. The scale p should be large 
enough in order to justify the calculation of Cd in QCD perturbation theory. 
In practice, using the standard methods of the Feynman-diagram technique, 
an explicit separation of distances is impossible in the quark-loop diagrams. 
One is forced to take into account both the soft parts of perturbative diagrams 
and the long-distance condensate effects simultaneously. This yields a certain 
amount of double counting, which is, fortunately, in many cases numerically 
insignificant, because the condensate contributions turn out to be much larger 
than the soft "tails" of perturbative diagrams. Moreover, if one does not go be- 
yond the two-loop diagrams in Figs. |b-|d and uses the leading-order Wilson 
coefficients Cd for d ^ 0, it is possible to rearrange the OPE including the soft 
(p 2 < fi 2 ) contributions of the perturbative diagrams in the definition of the 
condensates, ffhis "practical version" of OPE is discussed in more detail in a 
recent review a 

The list of the operators with low dimension entering Eq. (|2^) starts with 

3 = (31) 

and 

O^G^G^, (32) 

whose vacuum averages are known as the quark and gluon condensates, re- 
spectively. It is important that in QCD there are no colorless operators with 
lower dimensions, d — 1,2. 

The operators with d = 5, 6 are 

05=K,yG a ">, (33) 

Ot = (#VV>)(#» , (34) 
where r r s denote various combination of Lorentz and color matrices, and 

Oq = fabcG a ^ v G b J / 'G CCT/i . (35) 

The vacuum averages of the above operators are, correspondingly the quark- 
gluon, four-quark and three-gluon condensates. The condensates with d > 6 
usually play a minor role in the most QCD sum rule applications and we will 
not consider them in detail. 

One should mention that, in addition to rather mild effects of the quark 
scattering over the vacuum fields, there exist specific vacuum fluctuations at 
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short distances ~ 1/ y/Q^, which absorb the whole momentum of the exter- 
nal quark current. These effects, known as "direct instantons," |] violate the 
condensate expansion. For the vector currents considered here, short-distance 
nonperturbative effects may appear El only at very high dimensions (d > 10) 
and do not play any role in the truncated OPE. However, in the correlation 
functions of pseudoscalar (J p — - ) and scalar (J p = + ) quark and/or gluon 
currents, direct instantons are enhanced and important already at intermediate 
Q 2 - 

To proceed with the derivation of the QCD answer for 11^, on has to 
evaluate the Wilson coefficients of the condensate terms in (^0|). For illustra- 
tion, we demonstrate how the simplest contribution of the quark condensate 
is calculated. The relevant diagram is shown in Fig. 3d. One factorizes out 
all contributions to Tl^ containing one antiquark and one quark field, the 
remaining quark fields being contracted in the free-quark propagators: 

ng^fa) = * J d 4 xe^ x (0 | {r(x)^ (x, 0)7^(0) 

+^(0) lu S ji {0,x)^ i (x)} | 0) . (36) 

In the above, the fields tp and tp have to be treated as external vacuum fields 
with negligible momenta as compared with the momenta of the freely propa- 
gating off-shell quarks. In other words, it is possible to expand i>(x) and ip(x) 
around x = : 

if)(x) = <if>(0) + x" B P v(o) + .., 

#(x)=$(0)+$(Q)D P x p + .., (37) 

where D p is the covariant derivative, and the higher orders in this expansion 
are only relevant for the operators with d > 5. Substituting fl37|) in (^), one 
encounters the following vacuum matrix elements: 

(0 | | 0) = AS ij 5 al3 , (38) 

(0 I C D P % | 0) = B<F(7pW . (0 I C Dp ^(*) I 0) = BS»( 7p ) 0a , (39) 

where a, (3 are the bispinor indices, and the r.h.s. represent the most general 
decompositions in color and Dirac matrices, obeying color and spin conserva- 
tion. The constants A,B and B are easily obtained by multiplying both sides 

d They can be modeled employing the QCD instanton solution. Instanton-induced effects 
are discussed in the chapter by Shuryak in this book. 
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of Eqs. ( pq ) and ( |39j ) by S^Sa^ and 5 y (7 p ) Q( g, respectively, and taking traces. 
The result is : 

A=±{0\W\0), (40) 

| _ — ► 2777, — 

i?=-(0| V0^|O)=- — (0|^|0), B = -B. (41) 

In the last two relations the Dirac equation for the quark field $ ip(x) — 
—imil>(x) was applied. Substituting the expansion ( |37|) in ( |36| ) and using the 
expressions ( ^8| ) and (|3^) for the matrix elements, one obtains, after the inte- 
gration over x: 

n<T ) (?) = - «V) (0 I H I 0) , (42) 



yielding the Wilson coefficient 



C 3 (q 2 )= 2 -^. (43) 



The proportionality of the above expression to the quark mass is expected 
on general grounds. The quark condensate violates chiral symmetry and its 
contribution should vanish in the chiral limit m = 0. 

The derivation of higher dimensional terms of the OPE ( pp[ ) is more in- 
volved. A very useful tool, simplifying the calculational procedure, is the 
Fock-Schwinger gauge for the gluon field: 

{x - io)/'W = . (44) 

In this gauge, the gluon 4-potential A° is directly expressed in terms of the 
gluon field-strength tensor G*. This, and manjr other aspects of the cal- 
culational technique are explained in the review Ell serving as a very useful 
handbook for QCD sum rule practitioners (see also Ref. 15). 

The final result for the OPE, with all Wilson coefficients up to d = 6 taken 
into account, reads: 



n <« 2 > + 



4m 2 q 4 Uwq 4 

m 3 , — A a _„„ „ 2ira s r A a 7 A Q 



^{9s^u Y G a ^) + ((^5-^)^75-^)) 
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where the shorthand notation (O) = (0 | O | 0) is introduced. The above 
expression is valid for the light quarks ip = u,d, s. In this case the quark-gluon 
condensate contribution (diagram in Fig. ^|e) is suppressed by an extra factor 
m 2 /Q 2 . Note thatlhe three-gluon condensate term vanishes for the correlators 
of massless quarksE^I The accuracy of the OPE for II(q 2 ) is nnt limited by the 
above expression. Currently, n( per *) is known up to 0(a;?)[LH and the 0(a s ) 
corrections to many of the coefficients C n at d ^ are also available. If the 
quark is heavy (ip = c, fa), the quark condensate terms are suppressed | and one 
has, to d — 4 accuracyEl 

2 q 2 J Imn^( S ) (a s G a G a,iv ) 

^1) = -Jds s{s _ q2) + 4g?rg4 f(a), (46) 

4m 2 

where a = 1 — 4m 2 /g 2 , 

, 3(a + l)(a-l) 2 + l 3a 2 -2a + 3 
/W = 2^ ln V^T a* ' 

and Im# ert )(s) is given in pq). Note that at q 2 = the gluon condensate 
term has a suppression factor l/(4m 2 b ) 2 with respect to the perturbative part. 
The Wilson coefficients of d — 6, 8 terms in the OPE for the heavy quark 
correlator, corresponding to diagrams with three and four vacuum gLuons (e.g. 
the three-gluon condensate contribution) have also been calculated!!^ 



2. 7 What do we know about the vacuum condensates ? 

The vacuum condensates introduced in the OPE are purely nonperturbative 
parameters, hence, their numerical values (the condensate densities) cannot be 
directly calculated and have to be determined by other methods, to be briefly 
summarized now. 

The quark condensate plays a special role being responsible for the ob- 
served spontaneous breaking of the chiral symmetry in QCD (discussed in 
the chapters by Leutwyler, by Diakonov and Petrov, and by Meissner in this 
book). For this reason, the value of the quark condensate density was known 
long before it was used in QCD sum rules: 

f 2 m 2 

(W) = - 7T( V A - ~( 240 ± 10 MeV ) 3 ' ( 4? ) 
2(m M + m d ) 

e Heavy quarks do not develop their own vacuum condensates, being far off-shell at the 
momentum scale A vac , and the interactions of virtual heavy quarks with the light-quark 
condensates appear in higher order in a s . 
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for -0 = u,d and, in SU(3)-flavor approximation, also for ip = s. The above 
value corresponds to the normalization scale /i = 1 GeV.|] 

From first principles, very little is known about other condensates. Some 
attempts exist to calculate them on the lattice or in the models of the instanton 
vacuum. At present, it is still more safe, as it was done in the original workp 
to determine the condensate densities empirically, by fitting certain QCD sum 
rules to experimental data. Being universal, the condensates extracted from 
one sum rule can be used in many others. 

The gluon condensate density was originally derived from the sum rule 
( p0| ) for the correlation function of j M = c^^c currents. The hadronic spectral 
density was saturated by the experimentally known masses ray and decay con- 
stants fv of the charmonium levels V — J/ip,ip', employing quark-hadron 
duality (explained in the next subsection) for the heavier states. Substituting 
the correlation function ( [46| ) in the l.h.s. of (|20[), and fitting it to the r.h.s. at 
not very large n (in order to safely neglect d > 6 terms) the gluon condensate 
density, 

(— G" G aflv ) = (0.012 GeV 4 ) ± 30% , (48) 



has been obtained!] Note that the above quantity is scale-independent. The 
estimate ( fl8| ) has survived after many years, although claims urging to revise 
it ap pea r from time to time ip_the literature. In fact, an independent check of 



Eq. (48) has been carried outE3 by considering SVZ sum rules for two different 
correlators in the pion channel, one of them sensitive to the gluon condensate 
and the other one to the quark condensate. 

The quark-gluon condensate is "invisible" in the sum rules for vector 
mesons and has to be extracted from the correlation functions for baryon 
currents (they are considered in Sec. 3). The conventional parametrization 
is 

MvyG°' , ^)=^(#) ) (49) 

where the numerical value of ra\ has been estimated long ago and is still 
accepted: 

ml(l GeV) = 0.8 ± 0.2 GeV 2 . (50) 

Four-quark condensates with different combinations of r r matrices (one 
of the corresponding diagrams is shown in Fig. ||f) can be treated in the fac- 
torization approximation El which relies on the dominance of the intermediate 



f The condensate density is logarithmically dependent on the normalization scale if the 
underlying operator, in this case ipip, has a nonvanishing anomalous dimension. A natural 
choice is the scale /t separating the long and short distances in OPE. 
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vacuum state and allows to reduce each separate four-quark condensate to the 
square of the quark condensate. The general factorization formula reads: 



(^r r #r s V) = ^-^{(TrT r )(TrT s ) - Tr{Y r T s )}{^) 2 . (51) 

Applying it to the four-quark condensate terms in ( [45] ) one obtains, instead 
of the two terms in the square bracket, a more compact expression ^p(V'V') 2 • 
The product of the four-quark condensate and a s in (^5|) has a negligible scale- 
dependence. 

The three-gluon condensate contribution is not important for most of the 
interesting sum rules. An order of magnitude estimate of its density based on 
the instanton model istH 

(g^fabcG^G^G^) ~ 0.045 GeV 6 . (52) 

In the context of a recent development, let us mention the estimates of d = 7 
condensates obtained, in the model of the instanton vacuum Eil and using a 
factorization ansatztj 

One should admit that the accuracy of the condensate densities obtained 
from sum rules and/or invoking factorization is not very high and there is still 
room for a considerable improvement. Because of that, one cannot fully bene- 
fit from the improved Wilson coefficients available from perturbative calcula- 
tions. On the other hand, as we already noticed, perturbative loop diagrams 
include regions of soft momenta. Therefore, any update of condensates has to 
be combined with a more delicate procedure identifying and separating these 
regions in the perturbative coefficients, in spirit of the "practical version" of 
OPE. Regardless to their particular relevance for SVZ sum rules, the vacuum 
condensates are important nonperturbative characteristics of QCD in general, 
thus deserving dedicated studies. 



2.8 Use of quark-hadron duality 

Let us continue our derivation. Performing the Borel transformation of the 
Eq. (p5|), the QCD answer for II(M 2 ) in the sum rule ([l9]) can now be obtained. 
The result reads: 

oo oo 
f 2 e -m> v /M> + J ds p H {s)e -s/M> = _i_ L + «!&)\ J ds e -s/M* 





M 2 12M 2 81 M 4 
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(53) 



In the above, the four-quark condensates are factorized and a s is taken at the 
scale M | 

In addition, it is possible to estimate the integral over the excited and con- 
tinuum states in Eq. ( p^ ) using the following arguments. In the deep spacelike 
region q 2 — > — oo, where all power-suppressed condensate contributions can 
safely be neglected, the limit II(g 2 ) — > H^ pert ^ (q 2 ) is valid yielding an approxi- 
mate equation of the corresponding dispersion integrals: 



., / , ImII(s) 2 /j M^fs) 



i / ds - ( — ~k - q / ds —( — ' ( 54 ) 

s(s-q 2 ) J s{s-q 2 ) 



4/n- 



(again at q 2 — > — oo). In order to satisfy the above relation, known as the 
global quark-hadron duality^ the integrands on both sides of it should have the 
same asymptotics: 

Im n(s) -> Im n (pert) (s) at s -> +oo , (55) 

where, in our case, Im IT(s) is given by Eq. ( [Tl| ) and Im n( per *)(s) by Eq. (p8|). 
One could still allow Im II(s) to oscillate around the perturbative QCD limit. 
Eq. (|5f]) is an example of the local quark-hadron duality. From ( p4| ) and (|55|) it 
is postulated that at sufficiently large Q 2 = — q 2 the following approximation 
is valid: 

1 / ds 2^ ~ ~t / ds — ( 2V ' 56 

J s(s-q z ) n J s(s-q 2 ) 

where so is an effective threshold parameter which does not necessarily coincide 
with Sq. After the Borel transformation one obtains: 

- — ' 00 
dsp h {s)e- s ' M2 ~ i / d S Imn(" eri »( S )e- s / M2 . (57) 



The relation (J56J) and its Borel transformed version (J57]) represent the quark- 
hadron duality approximation used in SVZ sum rules to replace the integrals 
over excited and continuum states. The threshold parameter sq which, in 



s To do it more precisely, the Borel transformation has to be applied to the logarithmic 
dependence of a s on the scale q 2 . 

h For a more detailed discussion of quark-hadron duality see the chapter by Shifman in 
this book. 
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general, has to be fitted, is expected to be close to the mass squared of the 
first excited state of V. 

The assumption ( |56| ) is certainly weaker than local duality, because it 
involves integrals and, in particular, it is insensitive to oscillations of p h (s) 
around Im n( per *)(s). Moreover, at very large Q 2 , for the positive definite 
correlation functions, such as (Q), the validity of (^6|) is simply a mathematical 
consequence of the global duality. On the other hand, employing the duality 



approximation ( [56| ) in sum rules at intermediate Q 2 (or, equivalently, Eq. (57) 
at intermediate M 2 ) one also relies on the validity of the local duality ap- 
proximation d56|), starting from a certain finite sq. The latter assumption is 
not a strict consequence of the asymptotic freedom. It is therefore fair to 
call Eq. (p3q) "semilocal" duality. In this situation, the Borel transformation 



performed in the sum rule (53) is of a great importance, because, due to the ex- 
ponential suppression of the integral on the l.h.s., the sensitivity to the duality 
approximation (|57]) of this integral is not high. 

It is also important to recall that quark-hadron duality was confirmed in 
the channels accessible in e + e~ annihilation and in r lepton decays, where 
the spectral density p h (s) of excited and continuum states was measured and 
the hadronic dispersion integrals were compared with their perturbative QCD 
counterparts. For example, using the experimental data on J/ip,ip', .. — * 
decay widths, the quark-hadron duality for the charmonium channel was shown 
to hold to good accuracycJ 

Using the duality approximation ( |57j ) in Eq. (|5^), one can simply subtract 
the integral over p h (s) from the perturbative part on the r.h.s. This is the 
last step in the derivation procedure outlined and explained in the previous 
subsections. Our goal is achieved: the SVZ sum rule for the parameters of the 
ground-state hadron V can now be written explicitly. For definiteness, let as 
choose the case V = p. It corresponds to the correlation function (||), where 
the 1=1 combination of u and d quark currents is taken: 

= \(u^u~d^d) (58) 

with the decay constant defined as in Eq. (|Io|): 

(p°(p) | jP | 0) = 4™p^ P) * • ( 59 ) 



The resulting sum rule reads: 

Jp Uvr 2 V 
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(60) 



M 4 12 M 4 81 M 6 

Sq being the duality threshold for the p meson channel. 

2.9 SVZ sum rules at work 

Historically, Eq. @) was one of the first successful applications of the method. 
Recently, this sum rule was reanalyzed^ highlighting many interesting theoret- 
ical details. Here we will follow a more pragmatic procedure, trying to reveal 
the predictive power and the actual accuracy of the method. 

We start with discussing the choice of the Borel parameter. The sum 
rule ( |60| ) is not applicable at too small M 2 because the missing terms with 
higher-dimensional condensates and, potentially, also the short-distance non- 
perturbative effects, all proportional to large powers of 1/M 2 , may become 
too important to be neglected. Usually, the low limit on M 2 is adopted by 
demanding that in the truncated OPE the condensate term with the highest 
dimension remains a small fraction of the sum of all terms. This limit keeps 
the convergence of the condensate expansion under control and guarantees that 
one does not introduce a large error neglecting the higher-dimensional terms. 
At too large M 2 the quark-hadron duality approximation cannot be trusted. 
Therefore, one also has to choose an upper limit on M 2 , so that the exponen- 
tially suppressed contribution of the states above Sq remains a small part of 
the total dispersion integral. The value of Sq is not completely arbitrary, being 
correlated with the onset of excited states in the channel of the current if.. 
According to the experimental data, the resonance activity related to the first 
excited states in the p channel shows up at s ~ 1.5 2.0 GeV 2 , hence Sq in 
this vicinity can be expected. We have checked that in the sum rule (^0|), in 
the range 

0.5 < M 2 < 1.2 GeV 2 , (61) 

the d = 6 four-quark contribution is less than 10% and, simultaneously, the 
s > Sq part of the dispersion integral is less than 30 % of the total r.h.s.. 
One should emphasize that in certain sum rules the Borel window similar to 



Eq. (61) simply does not exist, that is, the lower limit of M 2 overshoots the 
upper one. The channels where the sum rules fail are usually plagued by "direct 
instantons" , and therefore the physical reason of the failure is understandable. 

After the range of M 2 is determined, we can fit the decay constant f p and 
the threshold parameter Sg from Eq. ( |60|) by demanding the maximal stability 
of f p within this range. The mass m p is fixed by its experimental value.j] The 

l The SVZ method has enough predictive potential to yield also the p meson mass m p by 



25 




Figure 6: p meson decay constant f p calculated from the SVZ sum rule (BOh as a function of 
the Borel parameter M 2 , fixing the threshold at Sg=1.7 GeV 2 (solid lincJTFor comparison, 
the results at sg=1.5 GeV 2 and s =2.0 GeV 2 are shown by the lower and upper dashed 
lines, respectively. The central experimental value is indicated by boxes, the triangles on the 
M 2 axis mark the allowed range of the Borel parameter. 



resulting numerical values of f p are shown in Fig. | as a function of M 2 , at 
the central values of the condensates (f47j) , ( 48 ) and at a s (lGeV) = 0.5. Within 



the interval (|6l|), the maximally stable f p corresponds to sq = 1.7 GeV 2 . 

In order to assess the predictive power of the sum rule ( |60| ) , it is necessary 
to estimate the theoretical uncertainties. They originate from the following 
sources: 

(a) dependence on the Borel parameter. The sum rule is an approximate 
relation, therefore it is not surprising that the final numerical result for the 
constant parameter f p changes with M 2 . As can be seen in Fig. ^, the variation 
of f p within the Borel window is small, 

j-Borel = 21() ^ 21? MeVj (g 2 ) 

and one can claim the success of the sum rule (|6(]). Indeed, a large instabil- 
ity with respect to the variation of AT 2 would indicate absence of important 
condensate contributions or may cast doubt over the reliability of the duality 
approximation. Apparently, all values of f p in the interval ( |62| ) can be equally 
trusted as sum rule predictions, therefore, the variation within the Borel win- 
dow has to be included in the total theoretical uncertainty. 

combining the sum rule ( |6o| ) with the same sum rule differentiated over 1/M 2 . We do not 
discuss this procedure here, concentrating on the determination of hadronic matrix elements 
such as } p . 
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(b) inaccurate knowledge of the condensate densities. Varying them within 
errors indicated in Eqs. ( p7| ) and ( |48| ) we obtain a small, ±1% change of the 
prediction (|62|) . 

(c) neglect of the d > 6 terms in OPE. We assume that the neglected 
condensate terms are altogether not larger than the d = 6 contribution. This 
yields less than ±5% uncertainty in f p . 

(d) limited accuracy of perturbative contributions. Varying the scale of 
a s and adding the known higher-order corrections to the perturbative part 
one again observes a small, about ±3%, variation of f p . Absence of 0(a s ) 
corrections to d ^0 Wilson coefficients has a negligible impact on the result. 

Adding the individual uncertainties (a) — (d) linearly, which is a rather 
conservative attitude, we get about ±10% in total and obtain the following 
interval for the SVZ sum rule prediction: 

f p = 213 ± 20 MeV , (63) 

where the form of writing is just for convenience and there is no preferable 
central value. The uncertainties (c),(d) may somewhat decrease after future 
theoretical work. Note that the variation of Sq should not be counted as 
an independent uncertainty, since this parameter is fitted from the sum rule 
together with f p . On the other hand, it is helpful to have several alternatives 
for the quark-hadron duality ansatz, including the radial excited or continuum 
states and increasing the effective threshold. In the p channel one can try 
different models, e.g., p + p'+continuum with higher s(J, including also the 
finite widths of these resonances. The results for f p are practically very close 
to what one obtains with the simplest SVZ construction: p + continuum. Our 
estimate ( |63|) is in a good agreement with the experimental number 

j-exp = 216 ± 5 MeV (64) 

obtained from the measured leptonic width r(p° — * e + e~) = 6.77 ± 0.32 
keV. 

Another illustration of the method is provided by the SVZ sum rule for 
the pion decay constant f^, defined as 

(T(p)l4 7r) |0) = -^/ OT , (65) 

where 

jjP = m^ 5 d (66) 
is the axial- vector current interpolating the pion. The sum rulecl obtained from 
the correlation function of two currents looks very similar to Eq. (|so|) : 



fl = M 2 



r l 

4^2 



(l _ (l + ^) 
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Figure 7: The decay constant f n calculated from the SVZ sum rule at so=0.7 GeV 2 (solid 
line). For comparison, the experimental value is indicated with boxes, the triangles at the 
M 2 axis mark the allowed range of the Borel parameter. 
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but, importantly, differs in the sign of the four-quark condensate term. This 
difference between the nonperturbative interactions of the currents jjf^ and 

(it) 

jn with the QCD vacuum, to a large extent explains why the hadrons in the 
vector and axial channels are not alike, i.e. why there is a single p resonance 
in the first channel and a combination of pion with a\ meson in the secondul 
The sum rule prediction with the estimated uncertainty, 

U = 127 ± 15 MeV, (68) 

is in excellent agreement with the experimental number f w = 132 MeV. Many 
other '^classical" examples of the SVZ method can be found in the original 
papers u or in the reviewsu 

There are two important messages from the above analysis. The first one 
is that QCD sum rules are (and will remain) approximate and their accuracy 
cannot be improved beyond certain limits. The second message is more en- 
couraging. Within the SVZ method one is able to estimate the theoretical 
uncertainty of the predicted hadronic parameter. Similar estimates are impos- 
sible, e.g., in many quark models of hadrons, where the inputs are nonuniversal 
and have no direct relation to QCD. 
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3 Applications and development of the method 

In the past two decades, QCD sum rules have been applied to many problems of 
hadron physics, with an accuracy, in several cases, substantially improved with 
respect to the original analyses. An incomplete list of applications includes: 

• determination of the light (u, d, s) and heavy (c, b) quark masses; 

• masses and decay constants of light and heavy mesons and baryons; 

• form factors of mesons and baryons ; 

• valence quark distributions and spin structure functions of the nucleon; 
structure functions of the photon, p meson and pion; 

• hadronic matrix elements relevant for the description of K°—K°, Bd—Bd, 
B s - B s mixing; 

• strong couplings and magnetic moments of mesons and baryons; 

• calculation of the parameters of effective theories, such as chiral perturba- 
tion theory (vPT), heavy quark effective theory (HQET), nonrelativistic 
QCD (NRQCD); 

• spectroscopy and properties of non qq hadrons (gluonia, hybrids); 

• hadrons in nuclear matter, properties of hadronic matter at high tem- 
perature and density. 

It is hardly possible to include in one review an exhaustive detailed pre- 
sentation of all important results obtained in twenty years of investigations. 
For the interested reader, a lot of important information, not covered here, is 
accumulated in the papers listed in the bibliography. In this Section we focus 
only on a part of the topics listed above, with an emphasis on the results which 
seem to be particularly important for current and future experimental studies^ 
Our intention is also to discuss the possibilities to improve these results and 
to outline the problems that still need to be solved. 

Before starting the presentation, let us again emphasize that QCD sum 
rules are approximate relations. The predictions obtained from their analyses 
are characterized by uncertainties which are estimated by varying the input 
parameters (a s , condensates, Borel parameters, etc.) in the allowed ranges. 
For a convenient comparison with the results obtained by other methods, we 

J Fpfl an earlier survey, see the Appendix on QCD sum rules in the "The BaBar physics 
book"Ej Some of the results collected there are updated here. 
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will quote the sum rule predictions as: central value ± uncertainty. The latter 
cannot be considered as a statistical error. It just indicates the range of varia- 
tion of the result from the (central value — uncertainty) to the (central value 
+ uncertainty). 



3.1 Light quark masses: m„, mj, m s 

A precise determination of the quark masses is a task of paramount importance 
for the Standard Model and its extensions. Chiral perturbation theory allows 
to determine the ratios of the light quark masses: 

rrj Til III 

— = 18.9 ±0.8, — = 0.553 ±0.043, -y- - = 24.4 ±1.5, (69) 

m d m d ^{m u +m d ) 

obtained at the pext-to-leading order from the measured masses of the pseu- 
doscalar mesonscj The determination of the individual quark masses has at- 
tracted a lot of attention in the QCD sum rule community The correlators 
of the divergence of light quark (u, d) axial-vector currents have been stud- 
ied to calculate m u + ma- In addition, correlators of the divergence of vector 
and axial- vector strange currents have been analyzed to determine m s , using 
experimental information on Kit and Kitt: resonances. 

Let us present in more detail this set of results. For further use, it is 
important to mention that the quark masses depend on the normalization 
scale \i through the renormalization group equation: 

M^«i<z(m) = -7(««)"^«(M)- (70) 

The scale dependence can be expressed as: 

m q {y) = m q R(n), (71) 
where rh q is a renormalization invariant parameter and 

R{ ^^^ m ^ + ^_^Y^ +0{a ^ (72) 

7i and Pi are the first coefficients of the mass anomalous dimension 7 and of 
the QCD (3 function, respectively; in the MS scheme, for three colors and Nf 
flavors, they read: 

70 = 2 , 71 = ^ - ^N f , (73) 
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and 

2 19 
0o = 11 - -N f , fa = 51 - —N f . (74) 

In order to determine m u + md one can study the two-point correlation 
function 

* 5 (<Z 2 ) = ijd i x e^<0|TO 5 (z)jt(0)}|0) , (75) 

j5 being the divergence of the axial current (|66|): j§ — d^uy^j^d). The 
correlator ( |75| ) is particularly sensitive to the u and d quark masses, since 
h = ( m u + m d )ui^ 5 d. j_H_. 

In the perturbative contribution to ([75] ) up to 4- loop corrections E3E3 are 
available. In the condensate expansion terms up to d = 6 has been taken into 
account. The hadronic spectral function is expressed in terms of the pion pole 
and of higher resonances, and is constrained to satisfy the behavior predicted 
by chiral perturbation theory at the (37r) threshold. 

In Table [t]we present a set of the results for m u + rfid, in the MS scheme 
at fj, = 1 GeV (see also Refs. 28-30). In the same Table we also include recent 
lattice QCD determinations, which are typically provided for the combination 
(friu+md) /2 at /i = 2 GeV; the corresponding values at \x = 1 GeV are obtained 
using Eqs. (^2|)-(Q). One observes that, although the uncertainty has not 



(m u +m d )(/x = l GeV) (MeV) 


Rcf. 


comments 


15.2 ±2.0 


DR8Z0 


LO 


15.6 ±3.4 


N89E3 


it 


12.8 ±2.5 


P98ii 


0(a 3 s ) 


9.5 ± 1.6 
12.6 ± 1.3 
11.6 ± 1.1 


G97EL 
BGLMQDB 
AOOB 


(quenched) lattice QCD 

a 

u 



Table 1: QCD sum rule results for m u + (/i = 1 GeV), compared to the recent lattice 
QCD determinations (renormalized to the same scale). 



substantially improved in the years, the most recent sum rule determination 
of m u + rfid predicts a value smaller than the previous ones. This is related to 
the account of higher order QCD corrections and to the improvement of the 
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hadronic spectral function. Notice that, being combined with the chiral ratios 
in Eq. (^3{]), the results in Table [I] allow one to individually determine m u , ma 
and m s . 

A similar analysis can be applied to determine the strange quark mass. 
One considers the two-point correlator of the divergences of the strangeness 
changing vector current js = (s'y^u) : 

*(q 2 ) = i [d A x e^(0\T{ JS (x)jl(0)}\0) , (76) 



which is sensitive to m s , since js — i{m s — m u )su. 

It is worth describing this calculation in some detail. In QCD, the sec- 
ond derivative of the correlation function (f76|), \f r (q 2 ) = (d 2 /(dq 2 ) 2 )^/(q 2 ), is 
obtained from the dispersion relation: 

I" Jo ( s - Q ) 

The hadronic spectral function p(s) = — Im ^f(s) can be obtained inserting a 

set of intermediate states with strangeness |5| = 1 , J p = + and / = i in the 
correlator (fT6|), starting from the two-particle states: \Kir), \Krf), \Krf), etc. 
The \Kir) contribution to p(s) can be written as 

9 = 32^2 1 \ d (s)\ 0{s - s+), (78) 

where s± = (Mjf ± M w ) 2 . The function d(s) is the scalar K — > 7r form factor; 
in the low-s region it admits a linear expansion, with the parameters fixed by 
the one-loop chiral perturbation theory. Above the Kir threshold, the function 
d(s), hence p Kv (s), can be reconstructed assuming the dominance of the scalar, 
|5| = 1, i^(1430) and #£(1950) resonancesSS This procedure, however, 
only partially uses the experimental information on the scalar Kir system, 
where a sizeable non-resonant component is observedcJ Another possibility 
for reconstructing p K7T (s) consists in using the Omnes formula to account for 
information on the scalar Kir system from the measured I = i Kir scattering 
phase shift. With this procedure, the noxmalization condition of p K7r (s) at the 
Kir threshold is automatically satisfied^ 

Due to the positivity of the spectral function p(s), the procedures de- 
scribed above yield either a lower bound on m s , if only the low hadronic states 
are taken into account, or a determination of m s , if the rest of contributions 
is approximated by duality. Concerning the perturbative contribution to the 
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sum rule, early calculations included 0(a?) corrections; currently the O(a^) 
expression is uscdEao Also the condensates up to d = 6 are included-in the 
OPE. The results of this set of determinations are collected in Table [J 43 ! The 
main uncertainty in the determination of m s is due to the procedure of recon- 
structing the spectral function p(s)c3 

One could also use the correlator of the divergences of the strangeness 
changing axial curren±_proportional to m s + m u . The corresponding sum rule 
has been investigatedjScJ but the result for m s has a larger uncertainty caused 
by an insufficient information on K 7T7T resonances needed to reconstruct the 
hadronic spectral function. 



m s {n = l GcV) (MeV) 


Ref. 




comments 


171 ±15 


CDPS9I 


Jzl 


0{aj.) +rcsonances 


189 ±32 


JM95E 




a 


203.5 ±20 


CPS97 




0(al) +resonances 


140 ± 20 


cdnpa; 




0(al) + Omnes 


160 ± 30 


J98E 




it 


159 ± 11 


M99!. 


i 


FESR + Omnes 


155 ±25 


DPS991 


4ti 


axial current div. 


200 ± 40 ± 30 


CKP98 


47 


t decays 


164 ± 33 


PP99 4 




.. 


176 ± 37 ± 13 


KKP00 


| 


.. 


158.6 ± 18.7 ± 16.3 ± 13.3 


KMOol 






184 ± 26 


E97pi 




(unquenched) lattice QCD 


124 ±21 


G97Ef 




(quenched) lattice QCD 


139 ±9 


A99p 






138 ±5 


G99E 




« 


146 ± 12 


BGLMQ1 


ii 


« 


128 ±5 


G00^ 






144 ± 5 


A00^ 




.. 


117±9 


A00E 




(unquenched) lattice QCD 



Table 2: QCD sum rule results for m s (n = 1 GeV) compared to recent lattice QCD deter- 
minations (renormalized to /J, = 1 GeV). 



In Table |2j we also collect the outcome of the analysis of strange and 
nonstrange hadronic r decays. In this case, m s is determined retaining in the 
OPE the SU(3)-flavor breaking contributions, which depend on the strange 
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quark mass. The spread of the results is mainly caused by differences in the 
treatment of the higher-order perturbative terms. The sum rule results can be 
compared with recent lattice QCD determinations of m s , also given in Table |[ 
both in quenched approximation and for dynamical fermions. 

The compilation of results presented above clearly demonstrates the unique 
ability of QCD sum rules to analytically determine the light-quark masses, 
employing various correlation functions and systematically including, order by 
order, the contributions of the perturbative series and of the nonperturbative 
condensates. The main limitation is caused by an insufficient knowledge of 
the hadronic spectral function in the scalar and pseudoscalar channels. In 
principle, this difficulty can be solved when detailed experimental data become 
available. 



3.2 Heavy quark masses: m c , rrib 

The study of the charmonium system was probably one of the first applica- 
tions of the QCD sum rule methodtl The c quark mass can be determined, 
if one considers the two-point correlation function of two cry^c currents (dis- 
cussed in Sec. 2), and uses for the hadronic spectral density the experimental 
data on the e + e~ cross section into charm-anticharm states, in particular the 
precise information on the masses and the electronic widths of the J p = 1 ~ 
charmonium levels (J/V', ip', ■ ■ ■)■ So far, the sum rule was analyzed including 
the perturbative QCD two-loop corrections. The nonperturbative contribu- 
tion mainly depends on the gluon condensate (a s G 2 ), while the contributions 
of other vacuum condensates are believed to be small numerically. 



m c (GeV) 


m c (m c ) (GeV) 


Rcf. 




comments 


1.46 ±0.05 
1.42 ±0.03 


1.23ig;g| ± 0.03 


DGP9A 
N94E2 




Borel 
Borel, moments 




1.525 ± 0.040 ± 0.125 
1.33 ±0.08 
1.20 ± 0.04 ± 0.11 ±0.2 


APE98 
FNAL9S 
NRQCD9 


6C 

9B 


lattice QCD 
lattice NRQCD 



Table 3: QCD sum rule determinations of the c quark mass and recent lattice QCD results. 



In Table |3j we present two results for the charmed quark pole mass m c p3 
together with the recent lattice QCD determinations. The difference between 
the two analyses is that_in the first oneE3 the pole mass is directly computed, 
while in the second onet3 the MS running quark mass m c is determined, and 
then related to the pole mass. The two results are in agreement with each 
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other, and with the estimate m c (p 2 = —m 2 ) = 1.26 GeV (having about 
±0.1 GeV uncertainty) obtained in the original SVZ analysis. It would be 
important to update the sum rule determination of the charm quark mass, 
including the 0(a 2 ) pertiixbative correction L3 and the three- and four-gluon 
condensate contributions Ej 

In recent years, there has been an impressive progress in the determina- 
tion of the 6-quark mass, from the analysis of the two-point correlation function 
of Iry^b currents. In this analysis, the hadronic spectral function is obtained 
taking into account six T(nS) resonances, their masses and leptonic decay 
constants being precisely measured in e + e~ annihilation. The current activ- 
ity aims at working with the highest possible moments of sum rules, where 
NRQCD is a good approximation. In this framework one is able to perform 
the Coulomb resummation, taking into account relativistic and radiative cor- 
rections order by order. Notice that the gluon condensate contributions are 
negligible in the heavy quarkonium system. 

In Table ^ we present a summary ofLxecent determinations!^ together 
with an average of the lattice QCD resultsE3 From this Table we see that the 
preferred interval for the sum rule results is m& ~ 4.8 ±0.1 GeV. 

A detailed discussion of NRQCD and its application to the b quark mass 
problem is beyopd our task, and the interest reader is addressed to more spe- 
cialized reviewsl£j We only notice that there exists an alternative approach, 
based on taking the first few moments of the standard SVZ sum rules. In this 
case, the 6-quark mass is determined in a purely relativistic way. The drawback 
is the sensitivity to the tail of the hadronic spectral function corresponding to 
the open beauty production at energies larger than the resonance masses. This 
spectral function can be handled by employing the duality approximation, com- 
bining it with the experimentally measured inclusive e + e~ n -^ bb cross section. 
This kind of analysis was employed in the early papers □ and still deserves 
attention, having-jp mind the possibility of including the recently calculated 
0(a 2 ) correction^ which will improve the accuracy of the perturbative part. 

The task of an accurate evaluation of the b quark mass is very important. 
As we shall see in the following, a precise determination of nib will critically 
reduce the uncertainty in various sum rules for B mesons. 

3.3 Heavy-light mesons 

The decay constants /d (s) and /b (s) of the heavy-light mesons -D( s ), -B( s ) are 
defined by the matrix element of the heavy flavored axial-vector current 

< o\qi^ 5 Q\H( P ) >= if HPfl , (79) 
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m b (GeV) 


m b (m b ) (GeV) 


Ref.,-, 


Method 


4.72 ±0.05 
4.62 ±0.02 
4.827 ±0.007 
4.84 ±0.08 

4.88 ±0.10 
4.80 ±0.06 


4.19 ±0.06 

4.20 ±0.10 
4.20 ±0.06 

4.25 ±0.08 


DP9? 
N94 
V95 
JP99 
MY9 C 
H99 
PP98 
BS99 


1' 

57 

>e 

n 
J 

7( 
71 




standard SVZ 

a 

NRQCD 




4.26 ± 0.11 


H00E3 


lattice QCD (average) 



Table 4: The b quark mass from SVZ and NRQCD sum rules and the lattice QCD average. 

where Q = c,b, H = -D( s ), -B( s ) and q = u, c?(s), in the same normalization as 
/„• denned in Eq. ( |65| ) . An accurate calculation of these decay constants is very 
important for heavy flavor phenomenology, and QCD. sum rules were among 
the first analytical methods to predict fg and /yjEroO In these sum rules, 
perturbative 0(a s ) corrections and d < 6 condensates were taken in account. 

We would like to emphasize two aspects. Firstly, the result for fg is very 
sensitive to the value of the 6— quark pole mass. The latter is extracted from 
the analysis of the T system discussed in the previous subsection. Lowering 
the b— quark mass by 100 MeV increases fg by 30 to 40 MeV. The heavy 
quark mass dependence is less significant in the determination of fg. Secondly, 
the resummation of leading logarithmic contributions in the sum rule for fg , 
obtained in the framework of HQET, revealed thal./x,. has to be taken at a low 
energy scale, fi ~ 1 GeV, rather then at /i ~ m^Ou In full QCD with finite 
b— quark mass, the optimal position of the scale, which has to be determined 
from higher orders of QCD perturbation theory, remains an unsolved problem. 
Current sum rule calculations use a somewhat intermediate scale, equal to the 
Borel parameter: /i 2 = M 2 ~ 3 — 5 GeV 2 . The shift from ra b to lower scales 
produces another significant increase of Jb- In the case of Jd this effect is 
again milder. 

In Table §, we collect several results and compare them with the recent 
lattice QCD determinations. Notice that the account of the abovementioned 
effects has produced larger values for Jb than in the original calculations. 

The most recent sum rule results presented in Table || correspond to the 
following intervals: 



f D = 180 ±30 MeV, fs = 170 ± 30 MeV , 
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(80) 



Jd (MeV) 


}b (MeVJ 


Jxei. 


189 ± 49 
173 ± 22 

200 ± 20 


158 ±25 
168 ± 18 

180 ± 30 


DP8Z0 
D93L3 

KRWY9900 


195 ± 10±^ 
216 =h 11±| 
220 ± 3±| 4 


161 ± 16±^ 
173 ± 13±1 4 
218±5±| 1 


UKQCD99E1 

APEOO0 
UKQCDOO0 



Table 5: D and B meson leptonic decay constants and fg from QCD sum rules and 

lattice QCD. 

where the rather large uncertainties suggest that there is still a room for an 
improvement, as we shall argue below. 

The ratios JdJ /d and fsj Jb are calculated by including 0(m s ), SU(3)- 
flavor symmetry breaking corrections in the sum rules. In these ratios, the 
dependence on the heavy quark mass, as well as the effects of the radiative 
0(a s ) corrections, are less significant. A compilation of various determinations 
is presented in Table ^ again, one can try to summarize the results asl23 

f D Jfn = 1-19 ±0.08, f Bs /f B = 1.16 ± 0.09 . (81) 

The sum rule predictions for fo 3 and fo have to be compared with the 
results of the experimental measurements: fjj s = 280 ± 19 ± 28 ± 34 MeV 
and Jd = SOOt^t!^ MeV, respectivelyEj The data are still too uncertain to 
challenge theory. 

The decay constants of the vector mesons D* and B* , defined by the 
matrix elements: 

(0\q^Q\H*)=m H *f H *e^ (82) 
(Q = c, b, H* = D*,B*) have also been updatedEI 

fo* = 270 ± 35 MeV , f B * = 195 ± 35 MeV . (83) 

The corresponding matrix elements for the orbitally excited heavy mesons have 
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f If 
JD B 1 JD 


-f if 

IBs/ JB 


Ket. 


1.21 ±0.06 

1.15 ±0.04 
1.17 ±0.03 ±0.03 


1.22 ±0.02 
1.09 ±0.03 
1.16 ±0.04 
1.20 ±0.04 ±0.03 


D930,-. 

BCNP94E3 
N94^L 
HL96E3 


1.15±0.04±°°03 
l.ll±0.01±o 1 

i.o9 ±o.oij:8;8! 


i.i6±o.o6i8:jg 
i.i4 ±0.0218:! 
i.ii ±o.oil8:8i 


UKQCD99E1 

APEOO0 
UKQCDOO0 



Table 6: The ratios /d //u and fBs/fs- QCD sum rules (the four upper lines) versus 
lattice QCD predictions (the three lower lines). 

been determined, both for finite heavy quark massil and in HQET as discussed 
below in Sec. 3.9. 

Is it possible to improve the determination of the heavy meson decay con- 
stants? The analysis of the heavy- light quark systems does not allow an accu- 
rate, independent determination of the heavy quark masses, the latter should 
be obtained from the heavy quarkonium channels. However, once the accuracy 
of the b and c quark masses is increased, this immediately reduces the uncer- 
tainty of the heavy meson decay constants. The issue of a s corrections and of 
their proper scale and rcsummation requires further investigations, at least in 
the case of fg. In this case, the 0(a s ) corrections in the sum rule seem to be 
unusually large, and one may speculate that anomalously large soft parts of 
perturbative diagrams should be somehow separated. Furthermore, the con- 
tribution of the d = 7 (GGqq) condensate in the same sum rule is proportional 
to the heavy quark mass and could be sizeable, hence, it should be taken into 
account. Finally, improvement of the duality approximation for the hadronic 
spectral density is possible when enough experimental information on radially 
excited B and D resonances will be available. 

3-4 B c meson 

The studies of the charmonium system, in particular the determination of the 
masses and decay rates of J/t/), rj c and \c, were among the first applications of 
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the SVZ method and provided the first physical predictions of this approach. 
It is worth reminding the successful prediction of the mass of r] c a The anal- 
ysis was then extended to the bottomonium system (T,r]b,Xb )•••)■ With this 
rich experience in analyzing heavy quarkonia, QCD sum rules can be used 
to investigate the B c meson, the state with open beauty and charm observed 



recently by the CDF CollaborationE3 From the point of view of quark-gluon 
interactions, B c is intermediate between the cc and bb systems, and it shares 
with the two heavy quarkonia common dynamical properties. For example, it 
is possible to consider the heavy quark and antiquark as nonrelativistic parti- 
cles, and describe the bound state, adding then the relativistic corrections. On 
the other hand, B c , being the lightest hadron with open beauty and charm, 
decays weakly. Therefore, it provides us with a rather unique possibility of 
investigating weak decay form factors in a quarkonium system. 

In the framework of QCD sum rules, B c is investigated using the inter- 
polating current j 5 = i 6756. The analysis of the two-point correlator is very 
similar to that discussed in Sec. 2, including the gluon condensate and the 
0{a s ) correction, and summing up the Coulomb part of this correction in the 



the B c leptonic constant 
Clearly, the accuracy of 



nonrelativistic approximation. Various calculations 
yield results in the range f Be = 300 - 420 MeVEll 
this determination can still bc improved. This would be an important outcome, 
since the purely leptonic mode B c — > tv can be used to access the CKM matrix 
element V c b- Concerning the semileptonic B c decays, such as B c — ► J I %l)(rh^)tv 
and B c — ► B^lv, they have been investigated by three-point sum rulesfBEl 
following the .method we shall illustrate below. In particular, in this framework 
it is possible l3 to derive the relations among the form factors determined by 
the heavy-quark spin symmetry, i.e. exploiting the decoupling of the spin of 
the heavy quarks in the infinite quark mass limit. The numerical determina- 
tion of the form factors is still hampered by the absence of 0(a s ) corrections, 
which are currently estimated in the nonrelativistic Coulomb approximation 
at the zero- recoil point lj Abundant B c production is expected at hadron col- 
liders, and a carefuLexperimental investigation of this system will be possible 
in the near futureEl B c will represent an interesting testing ground for the 
QCD sum rule approach, and therefore the refinement of the theoretical pre- 
dictions concerning this system should be in the working plans of the sum rule 
practitioners. 



3. 5 Ioffe currents and sum rules for baryons 

QCD sum rules for baryons sugg ested by Ioffe provide an important demon- 
stration of the universality of the method, generalizing it from the quark- 
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antiqiiaxk states to the three-quark states. To construct the correlation func- 
tionsEM!3 one needs a baryon current, that is, a composite operator having the 
same quantum numbers as a given baryon. For the proton, several possibilities 
were studied: 

J N (x) = e abc (u aT (x)C %l u b (x)h 5 rd c (x) (84) 

or 

J' N (x) = e abc {u aT {x)Ca^u b {x)) lb a^d c {x) , (85) 

where a,6, c are color indices and C is the charge conjugation matrix. Other 
possible quark currents involve derivatives. Some criteria have to be adopted 
to single out the optimal interpolating current for a particular baryon. The 
first criterion is to choose a current with a minimal number of derivatives, in 
order to deal with low-dimensional spectral densities and, consequently, to min- 
imize the contribution of the excited states. Furthermore, the currents should 
maximize the projection onto the considered baryon state. This can be done, 
for example, by considering linear combinations of the interpolating currents, 
with the coefficient suitably chosen in order to maximize the overlap. Finally, 
it is possible to choose currents in such a way that the two-point functions 
are dominated by the perturbative contribution, with the condensate terms 
producing a hierarchical set of corrections. This last requirement suggests to 
use the current (pi|), instead of fl85|), to interpolate the proton. Studying the 
two-point correlation function 

H( ? ) = i J e^(Q\T{J N (x)J N (0)}\0) = ni(? 2 )+ dU 2 (q 2 ) (86) 

and neglecting the contribution of the continuum and of higher dimensional 
condensates, an astonishingly simple expression for the nucleon mass can be 
obtainedEil 

m 3 N ~-2(2Tr) 2 (qq)( l x = lGeV), (87) 

in agreement with experiment. 

With a little effort, one defines the interpolating currents for the L = 
baryonic octet: 

(88) 

(d aT (x)C^s b (x)) l5 r u c (x)] , 

and for the L — decuplet: 

J*(x) = e abc {u aT (x)C^u b (x))u c (x), 



J^(x)=e abc (u al (x)C lfl u b (x)) l5 rs c (x) , 



J a (x) = -e abc (s al (x)C^s b (x)) l5 ^u c (x) , 
J A (x) = \/|e 6c 



(u al (x)C^s b (x)) 75 rd c (x) 
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2(u aT (x)C^s b (x))u c (x) + {u aT {x)C ltl u\x))s c {x)\ , 

2(s aT (x)C^u b (x))s c (x) + (s aT (x)C^s b (x))u c (x)\ , 
tf(x) = e abc ( S aT (x)C^s b (x))s c (x). (89) 

The 0(a s ) radiative corrections to th^ two-point functions of baryon cur- 
rents are known for three massless quarkspl and have been workecLout also for 
the case of baryons containing one heavy and two massless quarkso This opens 
up interesting perspectives for more precise determinations of the baryon prop- 
erties, both in the light quark and in the heavy quark sector. Let us remind 
the reader that the analysis of the two-point correlator of baryonic currents 
was of prime importance for determining some characteristics of the QCD vac- 
uum. As a matter of fact, from the fit of the masses of baryons belonging to 
the L — octet and decuplet, it was possible to determine the values of the 
mixed quark-gluon condensate ( [49| ) and of the SU (3)-flavor breaking parame- 
ter 7 = (ss)/{qq) - 1 ~ -0.2& 

The extension to the case of baryons containing charm and beauty quarks 
is also straightforward. In the heavy quark sector, open problems concern the 
spectra of baryons containing more than one heavy (c and b) quark. These 
baryons will be experimentally observed and investigated at hadron colliders. 
Two-point correlation functkms of heavy-baryon currents allow one to predict 
the spectrum of these statesE3 For an overview of sum rule applications to the 
baryonic problems, and for an analysis of the features of various three-quark 
interpolating currents one can consult the available reviewsE3 

An important technique for the studies of the baryon dynamics is the 
method of external fieldsfH in which two-point correlators in the external static 
fields (such as magnetic field) are introduced. Originally, many static prop- 
erties of nucleons, such as magnetic moments, or nucleon matrix elements of 
axial-vector currents and other operators have been successfully calculated. 
The limits of this review do not allow us to discuss this technique and its ap- 
plications in mare detail. A thorough presentation can be found, e.g., in the 
lecture by IoffeQ To demonstrate the variety of problems that can be solved 
by employing the external field method or the closely related approach of the 
soft-pion field, let us mention some recent applications. In the light baryon 
sector, attention has been recently paid to the properties of negative parity 
baryons, such as A^*(1535). In particular, the experimental analysis has shown 
a suppression of the strong N*Nir coupling and an enhancement of the N*Nt] 
coupling. Interestingly, estimates of these couplings using the two-point corm* 
lators in the external (soft) pion field are in agreement with this observationE3 



^abc 



g ^abc 
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Other strong pion-baryon couplings have been analyzed by the same methodt 



3. 6 Three-point correlation junctions: form factors and decay amplitudes 

The method of QCD sum rules can be generalized in order to calculate the 
hadronic matrix elements of electromagnetic and weak transitions. In this 
case one starts from three-point vacuum correlation functions and uses dou- 
ble dispersion relations. This approach has been extensively used, both for 
light and heavy hadrons. The applications include the pion electromagnetic 



form factor 



i£3 radiative charmonium decays such as J/ip 



D 



an 



semileptonic and flavor-changing neutral current (FpNC) transitions cHl 
and, more recently, the radiative decays <j> — > (77, ?/)7t 08 
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In order to discuss the advantages and the difficulties of three-point sum 



rules, let us outline, as an example, the calculation 
nctic form factor defined by the matrix element: 



10c 



of the pion electromag- 



where q = p' — p and j e ™ is the electromagnetic current 



(90) 



(91) 



The starting point is the correlator of with two pion-interpolating currents 



T^x(p,p') = (i) 2 d i xd i ye^'^y\0\T{ ] ^\x) ] r{QW ) {y)m 



p'pAP + P')xT(p 2 ,p /2 ,q 2 ) + 



(92) 



where the momenta p, p' and q flow through the axial-vector and the electro- 
magnetic currents, respectively. In Eq. ( |9^ ) we have only shown the relevant 
kinematical structure, the others being denoted by the ellipses. 

Inserting in Eq. (92) two complete sets of hadronic states with the quantum 
numbers of the pion one obtains the dispersion relation: 



T(p 2 , P ' 2 ,q 2 



&^ 2 ) +/ dsds > AfiiO 

(ml -p 2 ){ml -p 12 ) J Rl2 (s - p 2 )(s' - p' 2 ) 



Pi 



(P 2 )[ ds^ 2+W 2 )f ds^L, (93) 
Jr.2 s'~p' 2 J Rl s-P 1 



where the double dispersion integral receives contributions from the excited 
and continuum states located in the region R12 of the (s, s') plane. The terms 
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containing the polynomials Pi and Pi arise from subtractions in the dispersion 
relation, and, therefore, two independent Borel transformations in p 2 and p' 2 
are needed to get rid of them. The same transformations enhance the double 
pole term with respect to the double integral in Eq. ( |9^ ) . 

On the other hand, the amplitude ( |92|) can be computed by a short- 
distance expansion, in terms of perturbative and condensate contributions: 

T(p 2 ,p' 2 ,q 2 ) = Y,C d (p 2 ,p' 2 ,q 2 ,v)(O d (n)) . (94) 
d 

The expansion is valid for large spacelike external momenta: \p 2 | , \p' 2 | 3> 
Ag CI) ; the squared momentum transfer Q 2 = —q 2 is also kept large in order 
to stay far away from the hadronic thresholds in the ^-channel. In Fig. we 




(<0 (f) 



Figure 8: Contributions to the three-point correlator Q94J) : 0(a a = 0) perturbative term (a) 
and some nonperturbative corrections (b— f). 

depict the perturbative diagram and several contributions of the condensates. 
In the chiral limit, the contributions of the quark and quark-gluon condensates 
vanish. 
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To proceed further, one has to match the expression ( |94| ) with the hadronic 
representation ( p3| ) . Invoking quark- hadron duality, one approximates the dou- 
ble dispersion integral in the region R\2 with the perturbative one. After that, 
the double Borel transformation inp 2 ,p' 2 is performed, introducing two corre- 
sponding parameters M, M', so that the sum rule for F^^Q 2 ) reads: 

flFAQ 2 ) = 4/ d S d S 'ImC°( S , S ',Q 2 )e-ifr-^ 



£^(M 2 , 

d=3 



M' 2 ,Q 2 ^)(O d ^)) 



(95) 



The domain R±2(sq), characterized by the threshold sj, is determined by the 
duality approximation and by the procedure of subtracting the contributions 
of excited and continuum states. More specifically, taking M' — M and con- 



side: 
rule 



100 



the contribution of condensates up to d — 6, yields the following sum 



ds ds' po(s, s', Q 2 )e 



48ttM 2 



52tt 
81M 4 



a s (^) 2 (l 



2Q 2 , 
13M 2 ' 



(96) 



where 



p (s,s',Q 2 



3Q 4 



1 



- 16^ 2 A 7 / 2 [ 3A(a 
Q 2 ) 2 - 4ss' and a - 



! ){<j + 2Q 2 )-\ 2 -5Q 2 {a + Q 2 y\ , (97) 



A = (s + s' 

The numerical analysis of (^6|) uses the values of condensates given in 
Sec. 2 and the threshold sj inferred from the two-point sum rule (|67|). For 
intermediate values of the momentum transfer, Q 2 = 1 -j- 3 GeV 2 the appro- 
priate range of Borel parameter is 0.7 < M 2 < 1.7 GeV 2 , where the sensitivity 
to the duality approximation is low and where there is a hierarchy of power 
corrections, the criteria adopted in the analysis of two-point SVZ sum rules 
in Sec. 2. For larger values of Q 2 the contributions of the higher-dimensional 
condensates containing terms ~ Q 2 /M 2 overwhelm the contributions of the 
low-dimensional terms. Therefore, the three-point sum rule (|96| ) cannot be 
used to predict the large Q 2 behavior of F 7T (Q 2 ). The reason of the failure can 
be traced back to the truncated local condensate approximation, which is too 
crude to reproduce the physical mechanisms governing the exclusive j*tt — > tt 
transition.| We shall return to this point in Sec. 4. On the other hand, there 

k As one remedy solving this problem it was suggested to use (model-dependent) nonlocal 
condensates t£j 
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is no doubt that the sum rule (^6|) is reliable in the region Q 2 = 1 3 GeV 2 , 
where also the numerical result Q 2 F 7r (Q 2 ) ~ 0.3 GeV 2 is in agreement with 
the experimental data (the latter are presented in Fig. |ll|). 

The calculation of the heavy meson form factors follows the same strategy 
as outlined above. Let us begin with some useful definitions. The matrix 
elements governing the weak transitions of D meson to a pseudoscalar P = K,n 
and vector V — K* , p final state are: 



(P(p')\qj»c\D(p))=f+ P (q 2 )(p+p%+fn P (q 2 )(p-p%, (98) 



and 



oyDV( 2\ 

( V(p', A)|^(l - K)c\D(p) ) = ^S-± e fi e (v)* P0P > 

ni£> + my ^ 1 

i{m D +m v )A^{q^> + iA * V } q2) (eW*.p)(p + p\ 
p mo + ray 



+ ^[Af( ? V^(9 2 )](^'P)(p-P%, 09) 
where is the polarization vector of V, 

A? v (, 2 ) = ^±^A?v {q 2 ) - ^-^A^(q 2 ) (100) 
Amy Amy 

and ^lf y (0) =A° V (0). 

Each of the form factors introduced above can be studied by introducing 
an appropriate three-point correlation function. For example, the form factor 
fr>K can De computed from 

T^(p,p')=(if[d i xd i ye i V-*-™\0\T{j^ (101) 



where D and K mesons are interpolated by the corresponding currents j^' — 
i C75U and j^' — u-jx^s, respectively. In this case, due to the presence of 
a large scale m c , the accessible region of the squared momentum transfer q 2 
includes also small positive values: q 2 <C m 2 . There is another important 
difference with respect to the sum rule for the pion form factor: In the OPE 



for the correlator (101), the most important nonperturbative effects are due 
to the quark and quark-gluon condensates, and are described by the diagrams 
in Fig. ||b,e. The contributions of the gluon condensate and of the four-quark 
condensates are negligible. The subsequent steps include the use of quark- 
hadron duality and the double Borel transformation in p 2 ,p' 2 . The Borel 
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parameters M, M' are now kept different, since they correspond to heavy and 
light mass scales, respectively. The sum rule has the form: 

m 2 D f D f K , + , 2 n -^1—^. 1 f , , , T / / 2 \ — — 4r 
Jdk\1 ) e ^ = — I as dslmC DK {s,s ,q ,fi)e ~ ~ 



+ ^C^ K (M 2 I M /2 , g 2 , M )(O rf (/i))- (102) 



d=3 



The explicit expressions for the Wilson coefficients Cfj K can be found in the 
literature E£2It£3 The threshold parameters Sq and Sq determining the domain 
Ri2(sq,Sq ) can be inferred, together with /d and fx, from the study of the 
corresponding two-point correlation functions. 







/5p(0) 


4^(0) 


^(0) 


V m (0) 


D -> if 

D — > 7T 


exp 110 
exp 111 


0.60 ±0.15 
0.76 ±0.03 

0.50 ±0.15 


0.50 ±0.15 
0.58 ±0.03 

0.5 ±0.2 


0.60 ±0.15 
0.41 ± 0.06 

0.4 ±0.1 


1.1 ±0.25 
1.06 ±0.09 

1.0 ±0.2 



Table 7: Form factors of the weak D — » P, V transitions at q = 0. The experimental 
numbers are obtained assuming the nearest pole dominance for the form factors. 



In Table [7] we present a set of results obtained at q 2 
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104,105 



The 



q dependence of the form factors can be predicted, although the procedure 
is non-trivial since one should stay far from the thresholds appearing in the 
double dispersion relation. The q 2 dependence of A\ and A2 turns out to be 
quite mild; on the contrary, the dependence of / + and V is compatible with 
the nearest pole dominance. The predictions, within their uncertainties, are in 
agreement with experiment. Note that the form factor f^ K is important for an 
independent determination of the CKM parameter |V^ S | from the semileptonic 
D — ► Klvi decaysa 

In order to predict the form factors of B transitions to light mesons 
(71", K, p), one has to re place c. pL.6 and D — > B in the abov e sum r ules yielding 
results for semileptonic 
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10' 



!£3 The method 



E3E5 and FCNC B transitions!] 
can be easily generalized to different processes and final states, including or- 
bital excitations, and has provided interesting predictions. An example is the 
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determination Ei3E2jll23 of the form factors relevant for the FCNC B 
and B -> K*£ + l~ decays: 

( K*(p', X)\sa tlv q u bn\B{p) ) = i^^'W'^il? 2 ) 
- ml-,) - ■ q(p + P%}T 2 (q 2 ) 



'Ik 



:(p + p)M(q 2 ), 



K 



(103) 



where b R = 1/2(1 + 75 )6, q = p - ft . The result MED 7i(0) = 0.35 ± 0.05 
allowed to predict the ratio T(B -> *f*7)/T(& -> s 7 ) = 0.17 ± 0.05, which 
agrees with the experimental measurements.^ 

However, a close inspection of the general structure of three-point sum 
rules for heavy-light transitions reveals a difficulty which is manifest in the 
parametric mj dependence of the various terms of the short-distance expansion. 
An example is the form factor Ai in the B — > p matrix element. In the 
limit of large nib, one observes that the coefficients of the quark and quark- 
gluon condensates grow with mj faster than the coefficient of the perturbative 
contribution. This is another manifestation of the difficulty in approximating 
the OPE expansion by the first few terms, revealed in the sum rule for the pion 
form factor at large Q 2 . The physical origin of this difficulty can be traced 
back to the mechanisms of pmducing the light state in the heavy meson decay, 



for large heavy quark massliia Of course, this problem could be irrelevant for 



the actual value of the 6-quark mass, and for particular processes and final 
states: an example is the transition B — > it. However, in order to avoid the 
problem ab initio, an operator expansion on the light-cone can be used to 
describe heavy-light transitions, as explained in the next Section. 

Concerning b — > c transitions, _B .decays to charmed states (B — > D^*\ D**) 
were investigated in early papersStij for finite c and b quark masses. Since 
that time, the common attention has shifted towards approaches incorporating 
the heavy quark flavor and spin symmetry, with the development of appropriate 
sum rules in HQET (see Sect. 3.9). Our opinion is that three-point sum rules 
for 6 — > c transitions and finite quark masses represent a viable approach com- 
plementary _ta HQET, which allows one to control the accuracy of the heavy 
quark limitEj Actually, in the finite-mass sum rules the important contribu- 
tions of hard-gluon exchanges in the diagram in Fig. @a are not yet available. 



' Recent results for the exclusive B 
0.34) ■ 10~ 5 and B{B+ -» K*+-y) = 
the most recent measurements give: B{b 
B{b -» si) = (3.11 ± 0.35 ± 0.80 ± 0.72) • 



K*"f transitions are: B[BB- — > K* ^) 



(3.76+°;^ ± 0.28) 



-0.83 



10" 



( 4 -55t ;L± 



For the inclusive_d 



10" 



(3.15 ± 0.35 ± 0.36 ± 0.26) ■ 10" 



ii( 



;cay, 
and 
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With the current progress of methods for computing many-loop diagrams, it 
should become possible to calculate these two-loop three-point diagrams with 
different quark masses. 



3. 7 Hadron structure functions 

A standard application of QCD to the hadronic structure functions is to study 
their logarithmic dependence on the momentum transfer Q 2 using the per- 
turbative evolution equations. The initial conditions for these equations are 
usually parametrized and fitted to the experimental data at some intermediate 
value of Q 2 . A direct analytical calculation of structure functions remains a 
challenging task. It is therefore very important that QCD sum rules are in 
a position to solve this problem, albeit approximately and within— a limited 
range of the Bjorken variable. The method was suggested by IoffellHil and was 
originally applied to the nucleon structure functions, to obtain, in particular, 
the valence u and d quark-parton distributions in the nucleon at intermediate 
Q 2 . The idea is to consider the four-point correlator 



Tjfi, = -i / <Fx d^y d^z e lq - x e 



iq-x ip-(z-y) / 



(0\T{J N (y)ff(x)j±(0)J N (z)}\0), 



(104) 

j = dj^(l — 75)?/ and j+ = wy^,(l — 7s)<i being the weak quark currents, 
and J N the nucleon interpolating current (pi}). Inserting complete sets of 





(a) 



(b) 



Figure 9: Four-point correlator in Eq. (104) (a) and its perturbative part (b) 



states with the nucleon quantum numbers and picking up the ground state, 
one encounters the forward scattering amplitude: 



Tau = ~i / d 4 x e^(N\T{ff(x)j±(0)}\N) 



(105) 
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which is directly related to the nucleon structure function -^(^Bj, Q 2 ) 
sured in deep inelastic neutrino-nucleon scattering. The key observation 



mea- 



ls 



that the imaginary part of the correlator (104) 

Im = — [T^p 2 , q 2 , s + ie) - T^ip 2 , q 2 , s - ^e)] (106) 

(s = (p+q) 2 ) can be expanded at short distances in the range of large Euclidean 
momenta p and q : Q 2 = -q 2 ~ 10 GeV 2 and \p 2 \ ~ 1 GcV 2 < Q 2 . For 
example, considering the lowest order perturbative diagram in Fig. |^b and 
taking the imaginary part of this diagram in s, one notices that the quark 
virtualities in the t channel arc of the order of p 2 XBj, with xsj = Q 2 /%P' q 
being the Bjorken variable. Therefore, the quarks are far off-shell if xsj is not 
too small. 



The OPE for the correlation function (104) is obtained E23 in a standard 
way in terms of perturbative and condensate contributions. In this calculation 
one retains only the leading powers in the ratio \p 2 \/Q 2 accounting for the 
leading twist in the structure functions (parton distributions). The dispersion 
relation is taken in one variable p 2 , which is the virtuality of the baryon current. 
One then performs a Borel transformation in p 2 and separates the contribution 
of the lowest-lying nucleon state from the contributions of the other poles and of 
the continuum. This is a rather sensitive point of the procedure, because, even 
after the Borel transformation, there is still a contamination of the ground- 
state term by "parasitic" contributions in the dispersion relation.^] The range 
of the Bjorken variable, which is accessible to the method, does not include 
small xsj — 0, where OPE breaks down. The region xgj — 1 is also out of 
the reach of this method, simply because it is the resonance region in the s 
channel. 

With all such caveats, the valence quark distributions in the nucleon have 
been computed: u v (xujl in the region 0.2 < xbj < 0.7 and d v (xBj) in the 



region 0.3 < xsj < 0.5LLH3 The accuracy of the calculation is at the level of 50%; 
within this accuracy the result is in agreement with the parton distributions 
extracted from experimental data. The mere fact that an analytical expression 
containing quark and gluon condensates reproduces rather subtle dynamical 
features of the nucleon is impressive and promising. 

Estimates of the polarized structure functions g\{xBj) and 32 {xbj) ob- 



tained by the same methodt£j can be presented as average values in the (rather 



m Raee«tly, in calculating the valence quark distribution in the pion, this problem was 
solved F 23 l by considering a more complicated kinematical configuration of the 4-point corre- 
lator, with two different virtualities of interpolating currents. In this case double dispersion 
relations and double Borel transformation are applicable, eliminating all "parasitic" terms 
and improving the accuracy of the sum rule. 
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narrow) accessible ranges of the variable xsj 







gi(x Bj )(0.b < x Bj < 0.7) = 0.05 ±50% (107) 



g2{x B j)(0.h < x Bj < 0.8) = -0.05 ± 50% . (108) 

These estimates are in a reasonable agreement with the data in the same 
intervals of the Bjorken variable: 



gi (x Bj )(0.5 < x Bj < 0.7) = 0.08 ± 0.02 , 125 (109) 
gi (x Bj ){0A < x B j < 0.7) = 0.08 ± 0.02 ± 0.01 , 126 (110) 



and 



g 2 (x Bj )(0.5 < x Bj < 0.8) = -0.037± 0.020 ±0.003% , 127 (111) 



Having an approximate but reliable method to calculate the structure func- 
tions at intermediate x B j, it is interesting to apply it to a structure function 
which is completely unknown and not even directly measurable in deep in- 
elastic scattering. This is the case of the chirality-violating structure functim: 
h\{xBj) which characterizes the transverse spin dynamics in the nuclcon 



L2S 



The dominant, :u=quark contribution to the proton structure function h\ has 



been estimated employing an appropriate four-point correlator. The pre- 



diction, in the region 0.3 < xbj < 0.5, is: 

hi(x B j)(0.3 < x Bj < 0.5) = 0.5 ± 50% . (112) 

The photon structure function F2 1 {x B j, Q 2 ) can also be obtainedE^ start- 
ing from the correlator of four electromagnetic currents with two different vir- 
tualities q 2 and p 2 . The imaginary part of this correlator, which determines 
the structure function of the photon with a virtuality p 2 , is calculated in a 
form of a condensate expansion, and in the approximation of small p 2 /Q 2 . In 
order to obtain the structure function of the real photon, a dispersion rela- 
tion in p 2 is used, expressing the virtual photon structure function in terms 
of the integral over hadronic states and employing the quark-hadron duality. 
The analyticity of the dispersion relation in p 2 permits the extrapolation to 
p 2 = 0. This method allows one to separate the hard and soft (hadronic) parts 
in the photon structure function. The hard part corresponds to the pointlike 
quark-photon interaction, whereas the soft (hadronic) part receives contribu- 
tions from the interactions of the photon with quark-antiquark states at large 
distances. The duality threshold Sq serves as an effective boundary of this sep- 
aration. The result for the structure function F2~f(xBj,Q 2 ), predicted in the 
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range 0.2 < x B j < 0.7 and at Q 2 ~ a few GeV 2 E^l agrees with the experimen- 
tal data. This agreement is impressive, because all inputs in this calculation 
are fixed by two-point sum rules. Furthermore, in this, approach it is possible 
to calculate the twist 4 corrections to F2 1 {x B j, Q 2 )pj and to find the gluon 
distribution in the photon 
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The four-point correlator with two heavy-light quark currents j'5 = Q"f§q 
and two electromagnetic heavy quark currents — Qj^Q ( Q = c or b) inter- 
polates a deep inelastic scattering on a heavy meson (D or B). The advantage 
in this case is that the correlator has a well-defined forward scattering limit, 
since the QQ states in the t channel are far from t — 0. Calculating this 
correlator in terms of the perturbative part plus condensate expansion and 
performing the double Borel transformation, it is possible to calculate a few 
first moments of the heavy quark parton distribution in D,B mesons(i^3 ob- 
taining, in a certain approximation, also the fragmentation functions of heavy 
quarks into heavy mesons. The calculation was done for finite quark masses, 
it would be interesting to repeat it in HQET. 

Finally, as a new direction in studying the structure functions, which is, 
however, out of the scope of our main presentation, one should mention the 
analysis of parton distributions in terms of a coordinate-space variable called 
the Ioffe time. These distributions have been obtained for proton and photon, 
employing the four-point correlation functions and sum rule techniqueliri 

3. 8 Matrix elements of effective operators 

The description of neutral meson oscillations involves a set of important hadronic 
parameters. The first one is -Bff(/x), defined by the matrix element 

(K°\(sj»(l - l 5 )d)(s^(l - l5 )d)\K°) = 2(1 + j-){f K m K fB K ^) (113) 

and representing the deviation of the K° — K° mixing amplitude from the 
vacuum saturation approximation. A similar definition holds for Bd, s , and the 



matrix elements analogous to (112) are expressed in terms of BB d , s (fJ-)- The 



leading dependence on the renormalization scale n is: 

B k (jm) =B K (a s (n))^ , B BdtS (fjL) = B Bd , s (a s (fij)* , (114) 

Bi being renormalization invariant quantities. 

There are two different methods to determine Bk- The first one is based 
on the analysis of the two-point correlator of the AS — 2 four-quark operator 
0\s = 2 — (s7 M (l — 75)d)(s7^(l — with the hadronic spectral function 
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receiving a contribution from the KK intermediate state. The second method 
consists in the calculation of the three-point function of Oas=2 and of the 
interpolating currents for K and K . The_determination of Bk from the two- 
point sum rule yields Bk — 0.55 ± 0.091H3 while three-point sum rules give 
B K = 0.4 - 0.9L! 



L3( 



Similar two- and three-point sum rules involving the relevant operator 
Oab=2 = (&7 M (1 — 75)<i)(67' 1 (l — 75)^) are employed to obtain the parameter 
Bs d - The results are compatible with the vacuum saturation approximation, 
namely Bs d — 1.00±0.15t23 All these calculations suffer from the uncertainties 
related to the neglect of nonfactorizable a s corrections. Notice that an accurate 
knowledge of the ratio r = fs d BB d / f% s BB a is nowadays indispensable for 
the analysis of the CKM unitarity triangle using information on Bd and B s 
oscillations. In this ratio, several uncertainties (from &-quark mass, radiative 
corrections, thresholds, etc.) should cancel out. Therefore, a direct extraction 
of the parameter r from the ratios of sum rules deserves to be worked out. 



3.9 QCD sum rules in HQET 

QCD sum rules are frequently applied in the framework of the Heavy Quark 
Effective Theory (HQET). One considers the correlation functions of quark cur- 
rents, where the heavy quarks are represented by their effective fields h v (x), v 
being the heavy quark four-velocity^ An example is the calculation of the pa- 
rameter F related to the B meson leptonic decay constant fg by the equation: 



Ib = C(m b )F 



A 



O 



(115) 



where the coefficient C(rrib) can be computed in perturbation theory. To deter- 
mine F, the SVZ method can be applied to the two-point correlation function: 

HH =ijd A x e**<0|T{^(20jM(0)}|0) (116) 

where u> = 2v ■ k and Jm{x) — hv{x)i"f5q{x) is the interpolating current of 
the pseudoscalar heavy-light mesons in HQET. Due to the heavy quark spin 
symmetry, F can also be computed from the two-point correlation function 
of the vector currents jv{x) = h v {x){ln — v^)q(x) interpolating heavy-light 
1~ mesons. The procedure starts from writing down a dispersion relation for 

"A more detailed discussion of HQET with the relation of the effective fields h v (x) to the 
heavy quark fields Q(x) can be found in the chapter by De Fazio in this book. 
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(116) in the variable u>: 



p2 poo n h {v\ 

II(w) = —J h / dv K ' + subtractions, (117) 

2A — u> J,. v — lu 

isolating the ground state contribution from the integral over the excited states 
and the continuum. The parameter A = ms-ms represents the binding energy 
of the light degrees of freedom in the heavy meson. The dispersion relation 
( |117| ) is then matched with the QCD expression, obtained for negative ui in 
terms of a perturbative term and condensate contributions: 



n{u)=U pert {u) + yC d ^4, . (118) 

' — ' — LJ ) a 



Finally, the Borel transformation is performed and, invoking quark-hadron du- 
ality, the contribution of higher states and of the continuum are approximated 
by the perturbative contribution above the threshold s/j. 

Another example is the calculation of the Isgur-Wise function At the 
leading order in the l/mb, c expansion, this function parametrizes the semilep- 
tonic B — >■ matrix elements. The form factox.£ ran be obtained from the 
sum rule for the three-point correlation functionOt 



E(w,u/,y) = i d A x dV ex^ k - x ~ k '-*'\0\T{jl,(x)M0)3 M (x')}\0), (119) 



y = v ■ v' being the product of the initial and final meson four-velocities and 
= hy'j^hy. Similar analyses can be applied to the form factors T\/%{y), 
r 3/2(j/)j etc, of B transitions to orbitally excited charm mesons. 

There are several advantages in applying QCD sum rules in HQET. For 
example, the relatively simple form of the Feynman rules in the effective theory 
allows to compute two-loop radiative corrections, not only for two-point. rn r rp - 



lators, such as ( jllq ), but also for three-point ones similar to (119)li3t-!ilt-!-U As 
we mentioned in Sec. 3. 6, such corrections are not yet available for finite quark 
masses in the three-point functions. Furthermore, a systematic renormaliza- 
tion group improvement in the current correlators can be performed; this is 
important, for example, in determining the choice of the scale of a s corrections 
in the calculation of the HQET parameters. 

Let us briefly mention a few numerical results obtained in this frame- 
work_Ea' 
tionsBCi' 



±he,-p^rameters F and A defined in Eq. flllq ) the sum rule predic- 
are, typically, 



F = 0.40 ± 0.06 GeV* , A = 0.9 ± 0.2 GeV . (120) 
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The parameter A = 570 ± 70 MeVliia is consistent with the determination of 
rrifc from the bottomonium system.^ The analogous quantities A + and A' have 
been computed for the + ,l + and 1 + , 2+ -hiiavy-light doublets, respectively, 
with the results A+ ~ A' = 0.9 — 1.0 GeVEda Combining this result with the 
calculation of A one can infer a prediction of the mass of P-wave qQ states 
which is consistent with the experimental data. 

The mass splitting between the lowest- lying 1~ and 0~|-s.t 
to the chromomagnetic interaction parameter A2 ; the result EM 46 



in 



(0- 



0.46 ± 0.14 GeV 2 fits well to the experimental measurement. 
More involved is the determination of the HQET parameter Ai, related to 
the kinetic energy of the b— quark in the B meson, where the results still have 
a large uncertaintyi 146 li 147 l The role of nondiagonal contributions (i.e. matrix 
elements between different radial excitations) to the double dispersion relation 
has to be better understoodG-3 

The Isgur-Wise form factors £(?/) and ri(y) have been computed at the 



next-to leading order in a s Pdjliid the functions T3 (jj) and rs (y) are known at 
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together with tl 



e form 



1:;') 
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the leading order in the strong coupling constant 
factors of subleading terms in the inverse heavy quark mass expansion 

The analysis has been extended to baryons containing one heavy quark, 
with the determination of the binding energy Aa 6 , the kinetic energy of 
heavy quark in the heavy baryon and four-quark matrix elements on the A^ 
the Isgur-Wise form factor governing the semileptonic Af, — > A c trans itio n 
and the matrix elements of decays such as A;, — > plv and A& — > A£ + £ 



It would be interesting to enlarge the field of applications of HQET sum 
rules by analyzing also four-point correlators with heavy-light currents in order 
to determine the heavy and light quark distributions in the heavy hadrons. 



4 Light-Cone Sum Rules 

4-1 The basics of the method 



The method of light-cone sum rules (LCSR)II^3~Mj is a fruiifojLljibrid of the 
SVZ technique and the theory of hard exclusive processesli^Jli^a The basic 
idea is to expand the products of currents near the light-cone. This procedure 
involves a partial resummation of local operators and avoids certain irregular- 
ities of the truncated OPE in the three-point sum rules. In recent years, the 
LCSR approach proved very useful in calculating various hadronic transition 



°For a more detailed discussion of the definition of A and its relation to the b quark mass 
see the chapter by Uraltsev in this book. 
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form factors. Within this approach, one is able to take into account both hard 
scattering and soft (end-point)li^J contributions. 

While SVZ sum rules employ vacuum-to-vacuum correlation functions, 
the starting object of LCSR is different. One considers a correlation function 
which is a T-product of two quark currents sandwiched between vacuum and 
an on-shell state. E£3 The latter can be a light-quark hadron (pion, kaon, p, 
K* , nucleon) or a photon. A physical example is provided by the process 
e + e~ — * 7r°e + e~. The hadronic part of this reaction is a fusion of two virtual 
photons into a single 7r° via quark e.m. currents. The corresponding amplitude, 
shown schematically in Fig. [To|a, has the structure of a typical LCSR correlation 
function: 

F^(p,q) = i J d^xe-^i^ip) | T{j; m (x)jt m (0)} | 0) 

= e flva0 p a q F(Q 2 ,{p-q) 2 ), (121) 

where p is the pion momentum, q and (p — q) are the photon momenta, Q 2 = 
—q 2 , j e ™ is the quark electromagnetic current ( |9l| ) and F is the invariant 
amplitude encoding the dynamics of the process. The chiral limit is adopted 
and p 2 = m 2 = 0. 



To derive the LCSR, one has to calculate the correlation function (121) in 
QCD, in the region of large Q 2 and \(p — q) 2 \ and to use dispersion relation 
to match the result of this calculation with hadronic matrix elements. Let 
us explain this procedure in more detail. Using unitarity in the channel of 
the current jf, m with the momentum p — q, i.e., inserting the complete set 



of hadronic states, one derives a dispersion relation for (121) in the variable 
(p — q) 2 keeping the second variable Q 2 fixed: 



F^(p,q) = 2- 



(t°(p) I ir I p°(p - q)) (p°(p ~ q) I oT I o) 



m*-(p- qf 



1 fjImF^(Q 2 ,s) 



ds (122) 



In the above, the ground-state contribution of the p meson contains the hadronic 
matrix element determining the 7*/? — > 7r transition form factor multiplied by 
the p meson decay constant: (p°(p — q) | j" n \ 0) = (fp/V2)m p ei P ^*. The dis- 
persion integral includes the contributions of excited and continuum states at 
s > Sq . The coefficient 2 takes into account the contribution of ui meson which 
is approximately equal to that of p. Calculating the amplitude F^ in QCD, 
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p-q , , q 

(a) 



(b) 




Figure 10: Light-cone expansion of the correlation function (121) 



one then applies the standard SVZ technique: the Borel transformation in 
the variable (jp — q) 2 and quark-hadron duality. The resulting sum rule allows 
to obtain the 7*p — > 7r form factor and to determine its dependence on the 
momentum transfer Q 2 . 



The correlation function (121) can be calculated by expanding the T prod- 
uct of quark currents near the light-cone x 2 = 0. This expansion is different 
from the local OPE used before and, as we shall see below, incorporates sum- 
mation of infinite series of local operators. 

Before proceeding with the light-cone expansion, we first have to convince 
ourselves that at sufficiently large Q 2 = —q 2 and \(p — q) 2 \ the dominant part 



of the integrand in the correlation function (121) stems from the region near 



the light-cone x — 0. Importantly, the pion momentum p does not need to 



vanish. Hence, the invariant variable v = q ■ p 
large: 



(q 2 — (p — q) 2 )/2 can also be 



W\~\(p-q) 2 \ 



Q 2 » A| CD . 



(123) 



5G 



It is convenient to define the ratio £ 



£ = 2*VQ 2 



(124) 



so that the region (123) corresponds to finite values of this ratio, £ ~ 1. Con- 
sider a reference frame where the pion three-momentum p is finite but small 
as compared to the photon virtualities: \p\ ~ /x, \po\ ~ and (i 2 « Q 2 , za In 
this frame, go ~ Q 2 £/(4m) + O(a0 and the argument q ■ x of the exponential 



function in Eq. (121) can be approximated as: 



q-x = q x - q 3 x 3 



Q__£ 

4/x 



■x 



'Q 4 £ 2 
16/i 2 



Q 2 ^ 3 



4/y 



(x - x 3 ) —x 3 . 



£ 



In order to avoid strong oscillations of the integrand, one has to demand xq - 
x 3 ~ 4/x/(Q 2 £), and simultaneously x 3 ~ £/(2/i). These two conditions yield 

x 2 ~ (a* + 4 M /Q 2 £) 2 ~ x\ + 4/Q 2 + 0{^ 2 /Q A ) , 



and, hence, a; 2 ~ l/Q 2 — > in the region (123). At the same time, there is 
no short-distance dominance, because xq ~ x 3 ~ £/(2/z) 3> 1/ \/Cp, indicating 
that an expansion in local operators around a; = is not applicable]] 

To proceed, let us calculate the leading-order contribution to the light- 
cone OPE of the correlator (121) corresponding to the diagram in Fig. [To|a. 
For simplicity we consider only the u-quark part of the currents without the 
electromagnetic charge factor. Contracting the u quark fields in (121), using 
the propagator of the free massless quark 



iS o (x,0) = (0 | T{u(x)u(0)} | 0) = 
and transforming 7^7^ -i^aupYlS + 



i 

2ir 2 x i ' 
, we obtain 



(125) 



F^ip, q) = -ie^ ap / d 4 x^^e-^ x (TT°(p) I u(x)Yl5u(0) I 0) . (126) 



To investigate the structure of the nonlocal quark-antiquark operator in Eq. fll26| ), 
let us expand it in local operators around x = 0: 



u(x) lplb u(Q) = -x) r lpl5 d(0) 



(127) 



p One may still use the local OPE for the soft pion, p = (£ = 0). In this case, (121) 



reduces to a one- variable amplitude similar to (|2|) and the short-distance dominance follows 
from the arguments presented in Section 2.1. 
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The matrix elements of these operators have the following general decomposi- 
tion: 



(%°(p)\u D ai D a2 ••• D ar 7 P 7 5 w|0) =(-i) r p ai p a2 ---Pu r PpMr 

+ {-i) r g aia2 Pa 3 --P ar PpM' r + ■■■ ■ ( 128 ) 

In the above, the first term is totally symmetric and traceless (at p 2 = 0) and 
contains only 4- vectors. There are other terms containing one or more g aiak , 
one of them displayed explicitly. Substituting the decom posi tion (127) in (126), 
integrating over x and using the definitions (128) and (|124|) one obtains 



F(Q 2 ,(p~q) 2 ) 



r=0 



Q 



4 



— ' r(r 

-2 v 



1) 



-ML 



(129) 



Since the variable £ ~ 1 in a generic exclusive kinematics with p ^ 0, all terms 
should be kept in each series in this expression. We have explicitly confirmed 
the qualitative conclusion made above: the expansion of F in local operators 
cannot be truncated at any finite order. One has to take into account and to 
sum up an infinite series of matrix elements M r , M' rl ... of lo cal o perators. On 
the other hand, there is a distinct hierarchy on the r.h.s. of ( 129 ). The second 
term containing M' r and further similar terms are suppressed by powers of a 
small parameter 1/Q 2 as compared with the first term containing M r . A closer 
investigation reveals that the diffe renc e between the local operators entering 
the first and the second term in fll29| ) is in their twist. Twist is defined as 
the difference between the dimension and the spin of a traceless and totally 
symmetric local operator. The lowest twist of the operators entering (128) 
is equal to two, simply because the operator without derivatives has dimen- 
sion 3 and Lorentz spin 1. Furthermore, after taking the matrix elements, the 
twist 2 components of the operators contribute only to the first, symmetric 
and traceless term of Eq. (128), containing M r . After multiplying both parts 



of Eq. ( |128| ) by g ai a 2 it becomes clear that the matrix elements M' r receive 
their contributions from the twist 4 operators, the lowest-dimension operator 

being u{D) 2 J_q15U. We conclude that one has to treat the nonlocal opera- 

in components 

= 0), the matrix 



tor in Eq. (126) by expanding it near the light-cone x 2 — 
corresponding to different twists. 

In the leading order of this expansion, at x 2 = (and p 2 



element in Eq. (126) has the following parametrization: 



{ir°(p)\u(x) lfMl5 u(0)\0} x 2= = 



-ip^^k I due lup ' x Lp^{u^) 
V 2 Jo 



(130) 
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where the function cp v (u, fx) is the pion light-cone distribution amplitude of 
twist 2, normalized to unity: J Q ^(u, fx)dix = l.| The normalization scale fx 
emerges due to the logarithmic dependence on x 2 and reflects the light-cone 
separation between the quark and antiquark fields in the operator. At x = 0, 



Eq. (130) is simply reduced to the matrix element (p5\) defining the pion decay 



constant. Furthermore, expanding both sides of Eq. (|130|) and comparing the 



l.h.s. with the expansions (127) and (|128|) we find that the moments of ip^iu) 



are related to the matrix elements of local twist- 2 operators: 

M r = — i— p f du u r (p n (u, fx). (131) 
v 2 Jo 

The function ip v (u), multiplied by /„., is a universal nonperturbative object 
encoding the long-distance dynamics of the pion. Together with the corre- 
sponding higher- twist distribution amplitudes, ip n (u) plays a similar role as 
the vacuum condensates play in SVZ sum rules. 



Substituting the definition (130) in Eq. (126), integrating over x, restoring 
the electromagnetic charge factor and adding the d-quark part we obtain the 
correlation function in the twist 2 approximation: 



F^(Q 2 ,(p-q) 2 ) = ^ jdv^jx) 



3 J uQ 2 — u(p — q) 2 



(132) 



F^ 2 \Q 2 ,(p~q) 2 ) = ^. I du ip„(u, li)T(Q 2 , (p - q) 2 , u, fi) (133) 



where u = 1 — u. Note that this representation has the form of a convolution 

l 

3 ~ 

o 

of the hard scattering amplitude T with the distribution amplitude cp w . The 
scale fx plays the role of the factorization scale which separates the contri- 
butions at x 2 < 1/fx 2 entering the hard scattering amplitude from the long- 
distance effects at x 2 > 1/fx 2 parametrized by the distribution amplitude. At 
zcroth order in a s , the hard amplitude 

T*HQ\(p-<lf,u) = uQ2 _ 1 u(p _ q)2 (134) 



q The complete definition includes the path-ordered factor necessary for gauge invari- 
da x fl A all (ax)\ a /2}, which is unity in the light-cone gauge, x^A af * = 



employed here. 
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is /i-independent. 

At this point one has to emphasize that the definition ( |130| ) is -actaially al- 
most ten years older than the method of LCSR. It was introducedE£3~t£a in per- 
turbative QCD studies of exclusive processes with large momentum transfer. 
In this approach, ip v (u) describes the distribution in the fraction of the longi- 
tudinal pion momentum carried by a valence quark in the infinite-momentum 
frame. The convolution ( |133| ) determining the asymptotic limit of th^hard 
exclusive 7*7* — > 7r° amplitude has also been obtained at that time. E£3 The 
0(a s ) correction to T has been calculated E£il from the diagrams similar to the 
one shown in Fig. [ic|b. 

We are now in a position to obtain a sum rule from the dispersion relation 
(122) matching it with the result of the light-cone expansion. Defining the 
matrix element 

(*°(p) I iT I p°(p - <?)> = ^(^K'^W , (135) 

in terms of the transition form factor F p7T (Q 2 ), we obtain, to leading twist 2 
accuracy: 

00 1 
V2f p FP*(Q 2 ) + t ±ImF(Q 2 ,s) = V2U f du <p v {u) 

m 2 -(p-q) 2 J s-(p- q) 2 3 J uQ 2 - u(p - q) 2 ' 



Furthermore, representing Eq. (132) in a form of the dispersion integral 

^V.fr-.ftjL ^fel (13 7) 
tt J S-(p- q) 2 



with 

-lmF^ tw2 \Q 2 , s) = J du^(u)6(uQ 2 - us) , (138) 

7T 3 J 



we obtain the duality approximation for the contribution of excited and con- 
tinuum states: 



oc 



ds ±lmF(Q 2 ,s) = r llm F^ 2 \Q 2 ,s) = y/2f« f du <p v (u) 

S - (P~ q) 2 J S-(p— q) 2 3 J uQ 2 - u(p - q) 2 ' 

s f> s p 

where Uq = Q 2 /( s o + Q 2 )- The duality threshold parameter Sq can be taken 
from the SVZ sum rule (pOj). Using Eq. (139), one can simply subtract the 



GO 



integral on the l.h.s. of Eq. (136) from the r.h.s. Performing the Borel trans- 



formation r _a r e finally obtain the LCSR for the form factor of the 7*p 



transition^ 

F " (Q2) - kl^ M ^M:^ + ^)- (140) 

The light-cone distribution amplitude (p^ , which is the necessary input for this 
sum rule, will be discussed below. 



To improve the accuracy of Eq. (140), one has to include into the sum 



rule not only 0(a s ) perturbative QCD corrections to the hard amplitude, but 
also higher twist effects. Physically, the latter take into account both the 
transverse momentum of the quark-antiquark state and the contributions of 
higher Fock states in the pion wave function. There are several sources of 
higher twist corrections in the light-cone expansion. First, one encounters the 
twist 4 contributions expanding the matrix element ( |130[ ) at x 2 — beyond 
the leading order: 







-iPnU / 


due mp - x 


Jo 










p-x J 



1 due m P- x g 2 (u) 7 (141) 

where gi2{ u ) are the light-cone distibution amplitudes of twist 4. Furthermore, 
the expansion of the quark propagator near the light-cone yields the matrix el- 
ements of the quark-antiquark-gluon operators corresponding to three-particle 
distribution amplitudes (diagram in Fig. |lO|c) . There are also four-quark con- 
tributions stemming from the propagator expansion, some of them shown in 
Fig. |l(i|d. These effects, together with other four-quark diagrams (Fig. [To|e) 
give rise (in the twist 6 approximation) to contributions which can be factor- 
ized into a product of the quark condensate and a two-particle distribution 
amplitude. This effect can be potentially important at intermediate Q 2 and 
has to be investigated case by case. Nonfactorizable four-quark contributions 
and all other higher twist terms suppressed by high powers of Q 2 can safely 
be neglected.]] 



Replacing the currents in the correlation function (121) by quark cur- 
rents with different spin-parity and flavor, one is able to interpolate other 
hadronic transitions involving the pion. This replacement will only change the 

' Note that twist 3 contributions to the correlation function (121) are proportional to 



and vanish in the chiral limit. 
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hard scattering amplitude. One and the same set of distribution amplitudes 
will enter the light-cone expansion of the new correlation function. The idea 
of combining perturbatively calculable, process-dependent coefficients with a 
universal long-distance input works here in complete analogy with SVZ sum 
rules. 



4-2 Light-cone distribution amplitudes 



As we have seen in the case of ipn(u), distribution amplitudes arc defined 
through the vacuum-hadron matrix elements of nonlocal operators composed 
of a certain number of quark and gluon fields taken at light-like separations. 
These operators emerge in the light-cone OPE of the T-product of currents. 
The relevant technique including the quark propagator expansion near .the, 
light-cone and the extraction of various twist components is quite general 
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It shares many common features with the technique used nowadays to study 
forward and nonforward deep-inelastic amplitudes. The distribution amplitude 
can be expanded employing the conformal symmetry of massless QCD. The 
conformal spin (partial wave) decomposition allows to represent each distribu- 
tion amplitude as a sum of certain orthogonal polynomials in the variable u. 
The coefficients of these polynomials are multiplicatively renormalizable, and 
have growing anomalous dimensions, so that, at sufficiently large normaliza- 
tion scale /i, only the first few terms in this expansion are relevant. The part 
of the distribution amplitude, which does not receive logarithmic renormaliza- 
tion is called asymptotic. The discussion of many important aspects of this 
analysis and of several interesting results obtained in recent years, is beyond 
the scope of this review, and we refer the reader to the literatureli£j'E££l Here, 
we only present the most important conformal expansion of the leading twist 
2 distribution amplitude: 



(Ptv(u, fx) = 6uu 



1+ £ a n (fi)C 3 J 2 (u-u) 

n=2,4,.. 



(142) 



where Cn / ' 2 are the Gegenbauer polynomials (for a derivation, see, e.g., Ref. 167). 
The coefficients a n are multiplicatively renormalizable: 



a n (ii) = a„(/z ) 



a s (po) 



7n/ 00 



and 



7n = C F 



n+1 



(n+l)(n + 2) 
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\k=l 



(143) 



(144) 



are the anomalous dimcnsionsli£3 At /x — > oo, a„(/i) vanish, and the limit 



o„ = corresponds to the asymptotic distribution amplitude 



(145) 



The values of the nonasymptotic coefficients a n at -a_ certain intermedi- 

for the moments 



ate scale /iq can be estimated from two-point sum rules 
J u n <p T (u, fi)du at low n. This method is attractive because it employs nonper- 
turbative information encoded in quark and gluon condensates. However, in 
practice the two-point sum rule determination of a n is not very accurate. One 
can also determine or, at least, restrict the nonasymptotic coefficients from the 
light-cone sum rules for measured hadronic quantities. Some examples will be 
discussed below. Regarding other methods, lattice OCT) has an almost unex- 
plored potential to calculate distribution amplitudes 



i£3 Furthermore, 



evaluated in the instanton vacuum model, at a low normalization scale 



16S 



In- 



terestingly, this calculation has produced a distribution close to the asymptotic 



ast 3 and 4 distribution amplitudes for tic. pian have also been worked 
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Use of the QCD equations of motion tSMZil yields rigorous relations 



outti 

between three- and two-body distribution amplitudes of the same twist, con- 
siderably reducing the number of independent parameters. 

To apply LCSR to hadronic matrix elements involving other pseudoscalar 
mesons (K,rj), one has to assess the SU(3)-violation effects in the light-cone 
expansion. There are several sources of such effects. The virtual s-quark prop- 
agator creates twist 3 contributions proportional to rn%, which differ from the 
analogous O(rn^) effects in the case of u, d quarks; the ratio of nonperturbative 
parameters fx I is larger than unity by about 20 % ; finally, the asymmetry 
between strange and nonstrange quark momentum distributions in the kaon 
has to be taken into account b y in troducing nonzero odd coefficients ai, as, ... 
in the Gegenbauer expansion (142) for <pk(u). 

The light-cone analysis for the vectox-tnesons (p, u, K*, 4>) has been worked 



out, including the update^ the twist 2p3 and a detailed study of twist 3 and 
4 distribution amplitudes 



taking into account the finite meson-mass effects. 
The distribution amplitudes for other than pseudoscalar and vector mesons 
still need to be determined. It is already possible to apply LCSR to hadronic 
matrix elements involving nucleons, since t he n ucleoli higher twist distribution 
amplitudes have been classified and_s±udiedtZ3 Finally, studies of the two-pion 
light-cone distribution amplitudes 1113 extend the field of a pplications of the 
LCSR method to the exclusive processes with two pionsr 7£ 



To complete our survey, let us mention the photon light-cone distribution 
amplitudes. They are important for many applications of LCSR to exclusive 



G3 



processes, where the photon is emitted at large distances. The photon dis- 
tribution amplitudes. were used, e.g., in the early study of the weak radiative 



transition £ — > pyliril In particular, the leading twist-2 distribution amplitude 



of the photon has the following definition: 



(j\u(x)o- a pu(0)\0) = e u (uu) / du(p 7 (u)F a p(uq,x) , (146) 

Jo 

where F a p(q, x) = (ef3q a — e a qp)e lq ' x is the photon field strength tensor, and the 
asymptotic distribution amplitude is ipy(u) — 6xu(l — u). The nonperturbative 
parameter \ — — 4 GeV 2 (at fi = 1 GeV) normalizing this distribution in units 
of the quark condensate density, can be determined from two-point sum rules 
with experimentally known hadronic parts. The update of this parameter and 
of higher twist photon distribution amplitudes is desirable. 



4-3 LCSR for the pion form factor 

An example of the application of LCSR is the calculation of th^pion electro 



magnetic form factor defined in Eq. (Bffl) . The original sum rulelilil was recently 



improved E_3 by calculating the 0(a s ) perturbative contribution of twist 2 and 



the factorizable twist 6 corrections. The starting object is the vacuum-pion 



correlation function similar to (121), but containing a pion interpolating cur- 



rent ju instead of jf, m ■ Theptesulting LCSR, at zeroth order in a 8 and in the 



twist 2 approximation, reads£2 



F ir(Q ) = I duip n (u,^ u ) exp ( ~-^p J — * Qi Jdsse 1 , (147) 



where ^(0) = — </4(l), and = Q 2 / (sj + Q 2 ), Sq is the duality threshold in 
the pion channel, taken from Eq. (|67|). The factorization scale fi^ = uQ 2 + uM 2 
corresponds to the average quark virtuality in the correlation function. This 



sum rule has a regular behavior at large Q 2 : the 1/Q 4 dependence of Eq. (147) 
at Q 2 — > oo corresponds to the soft end-point mechanism, provided that the 
integration region shrinks to the point u = 1. 

At 0(a s ), one recovers the leading ~ l/Q 2 asymptotic behavior corre- 
sponding to the hard scattering mechanism. Including this contribution in the 
LCSR and retaining the first two terms of the sum rule expansion in powers 
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Figure 11: The light-cone sum rule predictions for the pion electromagnetic form factor using 
asymptotic (dashed), CZ (dotted) distribution amplitude and the fit to the data (solid). Q 2 

is in GeV 2 . 



of 1/Q 2 one obtainsELZZ 



2a s 
3nQ- 



1 S 2 

~M^f du ^M + ipUo) fdsse-^ 



O 



.(148) 



The 0(1/(3 ) term in (148) coincidp-^with the well known expression for the 



asymptotics of the pion form factor 
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FAQ 2 



9Q 2 



da 



(149) 



obtained by the convolution of two twist-2 distribution amplitudes f n (u) of 
the initial and final pion with the 0{a s ) quark hard-scattering kernel. This 
coincidence can be easily traced provided that the SVZ sum rule (^) for f„ 
yields in the leading order ds e _s/ H!. = 4-7T 2 / 2 , and that J^duip% s (u)/u = 3. 



In addition to the twist 4 terms LZ3 the factorizable twist 6 contributions 



determined by the. quark condensate density, have been calculated and turned 
out to be smalltZZI Adding twist 4,6 terms to the twist 2 (leading and 0(a 8 )) 
parts yields the LCSR prediction for F 7T (Q 2 ) shown in Fig. [Tl]liZZl for the 



65 



pion distribution amplitude (142) with-bap choices: 02 = (asymptotic) and 
a 2 (lGeV) = 2/3 (the CZ distribution E3). For the twist 4 distribution am- 
plitude, their asymptotic expressions are taken, which nrovides a sufficient 
accuracy. The fit of the LCSR to the experimental dataliZS (with only 0,2 7^ 0) 
yields 



a 2 (lGeV) = 0.12 ±0.07; 



-0.05 
-0.07 > 



(150) 



where the first (second) uncertainty is experimental (theoretical). Within er- 
rors, this determination is compatible with the asymptotic distribution ampli- 
tude. 

Similar to the pion elastic form factor F„ , the LCSR for other transition 
form factors involving the pioi) have been obtained, including 7*p — > tt (chosen 



as a study case in Sec. 4.1) 
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and 7*7r — > aili_J Transition form factors of the 
pion to other light mesons are still unexplored within the LCSR method. 

Finally, LCSR contributed to the study of the 7*7 — > tt transition. The 



corresponding form factor is simply equal to the amplitude Eq. (121) at zero 
virtuality of one of the photons: F~ f7T (Q 2 ) = F(Q 2 , 0). Since the real photon is 
a long-distance object, the form factor F 7Tr (Q 2 ) contains nonperturbative con- 
tributions which are beyond the light-cone expansion of two electromagnetic 
quark currents. The leading 1/Q 2 asymptotics of F 1 ' R (Q 2 ) is well known 11^ 
and given by Eq. (U32) at (p — q) 2 — ► 0. The calculation of the contributions 



suppressed by powers of 1/Q is a nontrivial .task, at least when using three- 
point QCD sum rules and short-distance OPEE^J Within LCSR approach, one 
possibility not yet exploited is to employ the photon distribution amplitudes. 
Another way to estimate the 7*7 — > ix form factor is to use the hadronic dis- 
persion relation fll36| ) where the resonance term is determined by the LCSR 
(|140p and the integral over higher states is estimated using duality!!^ One can 
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3d that it does not 
again has a better 
of the asymptotic 
decreases the lead- 
ing order (twist 2,4) result by 15-20% leaving some room for nonasymptotic 
coefficients. The fit to the datajiields a 2 = 0.12 ± 0.03 at fj, = 2.4 GeV (if all 
other coefficients are neglected)t 



analytically continue this relation to (p 
contain subtraction terms. The result s 
agreement with the experimental data 
pion distribution amplitude. The 0(a s ) correction t£i 



18: 



Summarizing, the LCSR studies indicate that the twist 2 distribution am- 
plitude of the pion is close to its asymptotic shape, already at intermediate 
scales. More precise data on various form factors involving pion are needed to 
increase the accuracy of this determination. 
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Figure 12: Form factor of the 7*7 — » 7r° transition calculated with the asymptotic (solid), 
CZ (02(^0) = 2/3; long dashed) and BF(a2(^to) = 2/3, 0,4(^0) = 0-4; short-dashed) distri- 
bution amplitude of the pion (where ^0 = 0.5 GeV) in comparison with the experimental 
data points. The dash-dotted line corresponds to the Brodsky-Lepage interpolation formula. 

4-4 Strong Couplings 



The correlation function (121) can also be used to calculate the strong punr 
coupling, defined by the hadronic matrix element 



(ttVM-^I/ZV - q)) = e^eMeW^ g Ufm . 



(151) 



The idea is to employ the analyticity of the amplitude F(q 2 , (p — q) r) in two 
independent invariant variables q 2 and (p — q) 2 , the squares of momenta flow- 
ing through the two channels of electromagnetic currents. Inserting in these 
channels the complete sets of hadronic states with the quantum numbers of w 
and p, one obtains a double dispersion relation: 



Ffivip, Q) 



(7r (p)a;(-g)|p (p-g))(0|i« m |a;(-g))(p (p- g )|ir|0) 



■/ 

-R 12 



(ml 



ds\ds2 



q 2 )(m 2 



(p - q) 2 ) 



pL(si,s 2 ) 



(S2 - q )(si - (p- q) 2 ) 



+ {oj <-> p} 



(152) 



where p^ v is the double spectral density of the excited and continuum states 
and R12 is the region occupied by these states in the (si, S2) -plane. Eq. (152) 
is very similar to the double dispersion relation ( |93| ) used in Sec. 3 for the three- 
point correlation function. The terms arising from subtractions are denoted by 
ellipses. The ground-state resonance contribution includes t wo (a pproximately 
equal) combinations of p and u>. Substituting the definition (151) in the above 
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relation and expressing the matrix elements of the electromagnetic currents 
through the decay constants f PtU yields the dispersion relation for the invariant 
amplitude: 

F{q\ (p qf) = ^ ^f*=L ^ + ... . (153) 



To obtain the sum rule, one has to match Eq. (|153j) to the light-cone expansion 
(132), in the region of large \q 2 \ and \(p — q) 2 \- The double dispersion integral 



over continuum and excited states can be replaced by its dual counterpart, 
obtained by integrating the double imaginary part of F(q 2 1 (p — q) 2 ) over a 
certain region in the (si, S2)-plane. The lower boundary is = and the 
upper boundary is characterized by a single threshold parameter Sq (again 
assuming that Sq ~ Sq). The shape of this region is not very important because 
the double imaginary part of the amplitude Im Sl Im S2 F{s\, S2) obtained from 
( |132| ) is concentrated on the diagonal. su.— .S2- A detailed discussion of this 



procedure is available in the literature c£-i 
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Furthermore, one applies the Borel transformations with respect to the 
variables q 2 and (p — q) 2 , introducing two independent Borel parameters Mi 
and Mi-, respectively. The following formula is useful for this derivation: 



where M 2 = M?M$/(M? + M|) and u = M\I{M\ + Mf). It is natural to 
take M 2 = JWf = 2M 2 , i.e. uq = 1/2, having in mind almost equal mass scales 
in the p and ui channels. After Borel transformation, all subtraction terms 
vanish and the sum rule converts into a simple relation: 

- ^fe m ?/ M2 Af 2 (l-e- s o7Af 2 )^ (1/ 2) + 0( as ) + higher twists, (155) 

where we used f u ~ /p/3, and the subtracted exponent is the result of the 
quark-hadron duality approximation The twist 4,6 terms not shown explicitly 
have also been calculated^J whereas the 0(a s ) correction is not yet avail- 
able. The interesting point is that the sum rule for the strong coupling is 
determined by the value of the light-cone distribution amplitude at the middle 
point u — 1/2. This value is particularly sensitive to the nonasymptotic coeffi- 



cients a,2n- Therefore, the sum rule (155) together with analogous relations, for 



other measured strong couplings (like the LCSR for the ttNN coupling LLSi-a) is 
useful for constraining the nonasymptotic parts of distribution amplitude. For 
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the asym ptotic tp„ the numerical result g^, jm — 12.5 ±2 GeV" 1 obtained from 
Eq. ( |155| ) agrees with the experimental value g^p^ = 14 ± 2 GeV _1 EJ 

The methpxLpf external fields, briefly mentioned in Sec. 3.5, also allows to 



calculate g wp7r li£3 In fact, this method corresponds to the soft pion limit of the 



light-cone expansionE£j In this limit one can apply a short-distance expansion 
in terms of local operators with increasing dimensions. From that one may 
argue that sum rules in the external field can still be used for determinations 
of the strong couplings if the pion is soft, say, p w ~ 100 MeV. On the other 
hand, the external field technique is not applicable in cases when the pion has a 
relatively large momentum, like, e.g., in the p — ► irir decay, whereas LCSR work 
for any pion momentum. Moreover, in the one- variable dispersion relation used 
in the external field approach, the contributions containing transitions between 
the ground and excited states remain unsuppressed after Borel transformation, 
thus, reducing the accuracy of the sum rule. 

4-5 Heavy-to-light form factors and couplings 

During the last years, LCSR have been extensively used to predict the form 
factors of various transitions of heavy (B,D) to light Att, K, p, K*,...) hadrons. 
Reviews describing these studies are availablel 185 !'! 187 ! Here we only briefly out- 
line the main points, and present some recent results. 

Let us consider, as an example, the calculation of the form factor /^(f? 2 ) 
(the definition is similar to Eq. (]98|)). The underlying correlation function in 
this case is 



F»(p,q) = i J dxe^ x (n(p)\T{u(x)^b(x),m b b(0)ij 5 d(0)}\0) 

= F{q\ (p + qf)p^ + F{q\ (p + qf)q^ , (156) 

where the weak current uj^b is combined with the current bij§d interpolating 
the B meson. The leading order in a s corresponds to the diagram in Fig. |i"c| a 
with the virtual b quark. In the twist 2 approximation, one obtains, for the 
relevant invariant amplitude, 



dU (pir{u,p b ) 

uq 2 — u(p + q) 2 



Comparison of this expression with the result (|l32j) for the correlator of the 
light-quark currents nicely demonstrates the universality of the method: both 
expressions are determined by one and the same distribution amplitude ip T . 



The hard scattering amplitudes are different, the one in Eq. (157) being gen 



erated by the heavy quark propagator. Also the factorization scale of ip v in 
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Eq. (157) should be adjusted to the characteristic virtualities of the heavy 
quark. 

The LCSR is obtained by matching the light-cone exnansion of F with the 
dispersion relation in the B channel. The result readsE£jt 



188 



tiM 2 ) = 



i 



\1- 2 r I ''" 
1 m b U / — ex P 

u 



BJB 

0(a s ) + higher twists 



q 2 (l-u) 



uAP 



(158) 



where Uq 



m| — q 2 )/(sQ — q 2 ), s B being the duality-threshold parameter 
in the B channel. This sum rule is valid at q 2 , sufficiently lower than m 2 , in 
order to stay far away from the hadronic states in the channel of the weak 
current. In contrast to the case of the LCSR for the pion form factor, the twist 



3 contributions, not shown explicitly in Eq. (158), are important. 

Importantly, LCSR (158) has a regular behavior in the m b — * oo limit. 
If one employs the scaling relations for mass parameters and decay constants: 



m B = m b - 



•A, 



' + 2m b u , M 2 = 2m b r, f B 



/tub, where A, luq, t 



and Jb are independent of m b , the LCSR can be expanded in the heavy mass. 
The higher-twist corrections either scale with the same power of m b as the 
leading-twist term, or they are suppressed by extra powers of m b . Furthermore, 
the asymptotic scaling behavior sharply differs at small and large momentum 



transfers: at q 2 = one has / + (0) 
X does not scale with m b , f + (p 2 ) 



-3/2 L5 



1/2 



and at p 



2m b x, where 



i23 Hence, the large p 2 behavior is 



consistent with the asymptotic dependence predicted in HQET for the region 
of small pion momentumE23 Another important feature is that, similar to case 
of the pion form factor, the heavy-to-light form factors calculated from LCSR 



receive cc 
transition 



nlributions from both soft and hard mechanisms of the exclusive 
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twist 2 term 



192 



193 



At the rnrrpnt accuracy, Eq. (158) includes the 0(a s ) correction to the 



am 



gluon contributions 



INN 



L8 1 



twist 3 and 4 quark-antiquark and quark- antiquark- 
Concerning the numerical analysis of this LCSR, one 
has to mention that the sensitivity of f + {q 2 ) to m b is considerably reduced if 
Jb and Sq are taken from the two-point sum rule considered in Sec. 3.3. There 
is also a spectacular cancelation between 0(a s ) corrections to these two sum 
rules. Furthermore, the sensitivity to the light-cone distribution amplitudes 
is also not high because in the sum rule (158) these normalized distributions 
are convoluted with relatively smooth coefficient functions over a wide region 
of u starting from ~ 0.5 to 1. In particular, the result obtained from the 



sum rule (158) is less sensitive to the values of nonasymptotic coefficients in 
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the distribution amplitude tp^, than the LCSR result for the pion form factor. 
A typical overall uncertainty of heavy-to light form factors obtained by this 
method amounts to 15-20 %, if one employs constraints on nonasymptotic 
coefficients a% 4 .. provided by the comparison of the LCSR for the pion form 
factors (F w , F 17r ) and for the strong couplings with experimental data. To 
somewhat reduce this uncertainty, one needs better determination of these 
coefficients. It is also important to calculate the 0(a s ) correction to the twist 



3 contribution in Eq. (|158|) 



The same correlation function ( |156[ ) can be used to estimate the B*Bn 
coupling, defined as (B* (q)-K~ (p) \ B~(p + q)) = —gB*Bw(p • € ' yB " %> ), from the 
double dispersion relation ip B and B* channels, as explained in Sec. 4.4. The 
corresponding LCSR readst^J 

1 m 2 + m 2 / m 2 b B \ 

JB'9b*Btt =—5 7—? — e ^ m b f 7T M [e^-e m* ^(1/2,^) 



+ 0(a s ) + higher twists . (159) 
Here again the NLO accuracy in twist 2 has recently been achieved El In 



Fig. 13L3 the form factor f^^{q 2 ), obtained using Eqs. (158) and ( |15 
displayed in comparison with the recent lattice results (see references 
original paper L3) . Other applications of LCSR include the B — > K 



11L. 


the 


188 


|l94 
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195. L96 



Fig. 



[||196 

* K* 



B — > p, K* , (j) form factors of weak and FCNC transitionsE 
shows, for example, the recent predictions on the form factors of B 
transition. It is easy to convert the sum rules for B mesons to the corre- 
sponding sum rules for D mesons by replacing b with c, and B^ with flW 



in Eqs. (158) and (|159| ). The D — > -k,K form factors (the latter being. . 
sensitive to the SU(3)-flavor violation) calculated within this methodE2J'E£SE3 
are in a good agreement with the experiment and with the lattice QCD. In 



9 


9D*Dtt 


9B'Btt 


method 


Ref. 


0.22 ± 0.06 


10.5 ±3.0 


22 ±7 


LCSR (NLO) 


KRWY990 


0.39 ±0.16 


9± 1 
11 ± 2 


32 ±6 
20 ±4 
28 ±6 


SP (LO) 

J? t 
11 


089 197 ._. 
CNDDFG95 F 98 l 
BBKR 184 


n 97+0.04+0.05 
u - z ' -0.02-0.02 






xPT+data 


S98 199 


0.42 ±0.04 ±0.08 






lattice QCD 


UKQCD98 m 



Table 8: Strong B* Bit and D*Dn couplings from LCSR and from the sum rules in the soft 
pion limit (SP) compared with some recent results of other methods. The scale-independent 
coupling is defined as g = firgH* HTr/Zm-H , H = B, D. 



Table || the sum rule predictions for the strong couplings qwd-k and gB*Bn 
are compared with the result obtained on the lattice, and with the recent pre- 
diction of xPT constrained by experimental, data. A compilation of results 

LCSR have been applied to ob- 
as well as the matrix elements 
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of many other approaches is also available 
tain the strong couplings <?d*_d p and gs*Bp 
for the radiative transitions between heavy mesons (B* — > Bj) E23 and heavy 



2111 



baryonstiiij The strong couplings of negative and positive parity heavy mesons 



and pions have also been determined, both in the heavy quark mass limit and 
for finite c and b quark masses T _with predictions on the widths of + , 1 + and 2 + 
charmed and beauty mesonsc°JlH23 LCSR with photon distribution amplitudes 
have been employed to estimate the long-distance contributions to B — > 
weak radiative decaysEH^ 

Finally it is important to mention the use of QCD sum rules for the anal- 
ysis of exclusive nonleptonic decays of heavy mesons. Nowadays, these decays 
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Figure 14: LCSR predictions for the B — > A"* transition form factors. The dashed lines 
indicate theoretical uncertainties, the points represent lattice QCD results. 
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attract a lot of attention, being one of the main objects for studying the mech- 
anism of CP-violation. The form factors and decay constants calculated from 
LCSR and SVZ sum rules can immediately be used for estimates of the non- 
leptonic decay amplitudes in the factorization approximation, in particular, 
for assessing the violation of the SU (3)-flavor (or [/-spin) symmetry in various 
B-decays. More subtle effects, such as violation of factorization, annihilation 
and penguins can, in principle, also be treated within the sum rule frame- 
work. First attempts to analyze nonfactorizable effects in B and D decays 
and n 
OPEE 

form factors of effective, operators emerging in the short-distance expansion of 
correlation functions E2j It would be extremely interesting to continue these 
studies trying to create a framework where the light-cone OPE is consistently 
employed at all stages of calculation. 



8E ,207120: 



d_kmg distance effects in B — > if* 7 decays employed short-distance 
ijjf208| xhe LCSR have only been partially used in order to obtain the 



5 Summary 

In this review we have tried to present a concise and updated picture of QCD 
sum rules and their applications. This was not an easy task, because after more 
than twenty years of development, the manifold of works using this method 
resembles a large tree with many branches penetrating into different fields of 
QCD and hadron phenomenology. Some of the branches are quite distant with 
respect to each other, and in order to find the relevant papers one has to search 
not only in the hep-ph, but also in the hep-th and nucl-th electronic archives. 
Nevertheless, the procedure used in all these applications is essentially one and 
the same, formulated in the original workEJ A) construct a correlation function, 
B) calculate it at some Euclidean scale using the operator product expansion 
and C) match it with the hadronic dispersion integral. Above, in Sec. 2, we 
have described this procedure step by step, deriving the SVZ sum rule for the 
decay constant f p as a study case. 

The applications of QCD sum rules considered in this review cover only a 
part of the work done in this field. Following our particular line of discussion, 
we could have missed some important references. The aim of our presentation 
was not just to demonstrate how the problems of hadron phenomenology are 
treated within this method, but to help the reader to assess the current status 
of QCD sum rules in general. Let us try to formulate the main conclusions of 
our discussion: 

— QCD sum rules, together with experimental data on hadronic spectral 
densities can successfully be used to determine quark masses and universal 
nonperturbative parameters, such as vacuum condensates; 
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—Within the sum rule framework, operating with a handful of inputs, 
one is able to reproduce many hadronic observables, such as decay constants, 
form factors and parton distributions, in a reasonable agreement with the 
experimental data; 

—The accuracy of the method is essentially limited. Nevertheless, since the 
correlation functions are field-theoretic objects defined in QCD and related to 
hadrons via rigorous dispersion relations, the theoretical uncertainties can be 
traced and estimated; 

—Today QCD sum rules are not limited by the "classical" SVZ approach 
based on two-point correlators. As we have seen, combining the sum rule 
technique with the light-cone expansion, it is possible to develop the LCSR 
method which avoids certain problems of the local condensate expansion and 
adequately describes QCD mechanisms of exclusive hadronic transitions. 

Concerning the open problems, the major one is to reduce theoretical un- 
certainties in the inputs and in the procedure. Here QCD sum rules alone are 
not sufficient. Interaction with other nonperturbative methods is very impor- 
tant, e.g. lattice QCD determination of condensates or light-cone distributions. 
Some new experimental results of hadron physics could be very helpful. As 
we have seen, accurate measurements of the form factors of light hadrons will 
allow us to fix the nonasymptotic coefficients in the light-cone distribution am- 
plitudes. Furthermore, a better knowledge of excited hadronic resonances with 
different quantum numbers could reduce systematic uncertainties introduced 
by the quark-hadron duality approximation. One may even suggest dedicated 
experimental studies. For example, semileptonic decays of charmed mesons 
can provide important information about the resonances in the Kir and Kitty 
system with different spin-parities. 

To outline the future perspectives of the method, let us give only one 
important example. In order to extract fundamental parameters of the elec- 
troweak theory and to search for new physics employing the current and future 
data on exclusive B and D decays, accurate predictions on the heavy meson 
decay amplitudes arc needed. Approximate QCD methods are the only tools 
we have at our disposal, and the analytical method of QCD sum rules has, 
undoubtfully, a large unexplored potential in this field. 
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